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RESUMO 


ALMEIDA, Renan Augusto Lisboa, Universidade Federal de Viqosa, Fevereiro, 2015 

KARDAR-PARISI-ZHANG UNIVERSALITY, ANOMALOUS SCALING AND 
CROSSOVER EFFECTS IN THE GROWTH OF CdTe THIN FILMS. Orienta- 
dor: Sukarno Olavo Ferreira. Co-Orientador: Tiago Jose de Oliveira. 


Neste trabalho estuda-se a dinamica de crescimento de lilmes finos de Telu- reto de Cad- 
mio (CdTe) para temperaturas de dcposigao (T) entre 150 °C e 300 °C. Uma relaqao entre a 
cvolugao dos morros e flutuagbcs de longos comprimentos de onda na surperffcie do lilme de 
CdTe e estabelecida. Encontra-se que escalas de curtos comprimentos de onda sao ditadas por 
uma compcticao entre o resultado da formacao de defeitos na borda de graos vizinhos colididos 
e entre um processo de relaxaqao originado da difusao e da dcposicao de partfculas (moleculas 
de CdTe) sobre essas regioes. Um modelo de Monte Carlo Cinetico corrobora as explicates. A 
medida que T e elevada, essa compcticao da origem a diferentes cenarios na escala de rugosi- 
dade tais como: crescimento descoiTelacionado, crossover de descorrelacionado para cresci¬ 
mento correlacionado e escala anomala transiente. Em particular, para T = 250 °C, mostra-se 
que flutuaqoes na superffcie de CdTe sao descritas pela celebre equaqao Kardar-Parisi-Zhang 
(KPZ), ao mesmo tempo que, a universalidade das distribuiqrocs de altura, rugosidade local e al- 
tura maxima para a classe KPZ e, finalmente, experimentalmente demonstrada. A dinamica das 
flutuaqoes na superffcie de filmes crescidos a outras temperaturas ainda e descrita pela equaqao 
KPZ, mas com diferentes valores para a tensao superficial (u) e para o excesso de velocidade 
(A). A saber, para T = 150 °C encontra-se um crescimento Poissoniano que indica v = A = 0. 
Para T = 200 °C, entretanto, um crossover aleatorio-para-KPZ e encontrado, com A > 0 neste 
segundo regime. A origem da escala KPZ para filmes crescidos a T e [200, 250] °C decorre 
da complexa dinamica de empacotamento dos graos durante a qual espacos nas vizinhanqas dos 
mesmos nao sao totalmente preenchidos. Esse mecanismo de agregaqao tern o mesmo efeito 
da agregaqao lateral do modelo de deposicao balfstica o qual leva a um excesso de velocidade 
(A > 0). Finalmente, para filmes crescidos a T = 300 °C demonstra-se que A < 0. Em partic¬ 
ular, o mecanismo KPZ para filmes crescidos a esta temperatura decorre da alta taxa de recusa 
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da deposigao de pa it feu las, a qual e dependente das inclinagoes locais. Este fenomeno pode ser 
explicado em termos do coeficiente de sticking , o qual e tao pequeno quao mais localmente in- 
clinada for a superffeie. Devido a efeitos de tempo-finito (crossover temporals e escala anomala) 
ocoiTendo em T = 200 °C e 300 °C, expoentes de escala falham em revelar a Classe de Univer- 
salidade do crescimento. Contudo, um novo metodo desenvolvido, que avanga sobre a simples 
comparagao entre expoentes e seus valores teoricamente esperados, permite-nos concluir que o 
crescimento de CdTe, em uma ampla faixa de temperatura, pertence a classe KPZ. 
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ALMEIDA, Renan Augusto Lisboa, Universidade Federal de Vi^osa, February, 2015 

KARDAR-PARISI-ZHANG UNIVERSALITY, ANOMALOUS SCALING AND 
CROSSOVER EFFECTS IN THE GROWTH OF CdTe THIN FILMS. Supervi¬ 
sor: Sukarno Olavo Ferreira. Co-Supervisor: Tiago Jose de Oliveira. 


In this work one reports on the growth dynamic of CdTe thin films for deposition temper¬ 
atures (T) in the range of 150 °C to 300 °C. A relation between the mound evolution and large- 
wavelength fluctuations at CdTe surface has been established. One finds that short-length scales 
arc dictated by an interplay between the effects of the formation of defects at colided boundaries 
of neighboring grains and a relaxation process which stems from the diffusion and deposition of 
particles (CdTe molecules) forward these regions. A Kinetic Monte Carlo model corroborates 
these reasonings. As T is increased, that competition gives rise to different scenarios in the 
roughening scaling such as: uncorrelated growth, a crossover from random to correlated growth 
and transient anomalous scaling. In particular, for T = 250 °C, one shows that fluctuations of 
CdTe surface are described by the celebrated Kardar-Parisi-Zhang (KPZ) equation, in the mean¬ 
time that, the universality of height, local roughness and maximal height distributions for the 
KPZ class is, finally, experimentally demonstrated. The dynamic of fluctuations at the CdTe 
surface for other temperatures still is described by the KPZ equation, but with different values 
for the superficial tension ( v ) and excess of velocity (A). Namely, for T = 150 °C one finds 
a Poissonian growth that indicates u = A = 0. For T = 200 °C, however, a Random-to-KPZ 
crossover is found, with A > 0 in the second regime. The origin of the KPZ scaling for films 
grown at T G [200, 250] °C stems from a complex dynamic of grain packing during which all 
available space at the neighborhood of grains are not filled. This aggregation mechanism has the 
same effect of the lateral aggregation of the balistic deposition model which leads to an excess 
of velocity (A > 0). Finally, for films grown at T = 300 °C one demonstrates that a KPZ growth 
with A < 0 takes place. In particular, the KPZ mechanism at this T comes from the high refuse 
rate of the deposition of particles, which depends on the local slopes. This phenomenon can be 
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explained in terms of the sticking coefficient which is so smaller as more locally inclinated is 
the surface. Due to finite-time effects (temporal crossover and anomalous scaling) taking place 
in T = 200 °C and T = 300 °C, scale exponents fail in reveal the Universality Class of the 
growth. Notwhithstanding, a new scheme developed, which advances over the simple compar¬ 
ison between exponents and their theoretically predicted values, allow us, surely, to conclude 
that the growth of CdTe, in a wide range of deposition temperature, belongs to the KPZ class. 



1 


Introduction 


1.1 Motivation 

Semiconductor thin films are the basis of our opto-electronic technology and can be found every¬ 
where [1]. Great part of the current thin-film stage was due the vacuum technology adventum 
at 1950’s and 1960’s, which supported the developing of sophisticated growth techniques as 
Molecular Beam Epitaxy (MBE), Hot Wall Epitaxy (HWE), Chemical Vapor Deposition (CVD) 
and others [1-3], Since then, the quality and control on doping, thickness, chemistry compo¬ 
sition and structure of films became so accurate that the results have been touching directly 
our life: computers each time smaller and more powerful (including cell phones and others 
mobile devices), systems of global position localization (GPS) available to population, lasers of 
several wavelengths, medical equipments of high social impact as magnetic resonance, positron- 
electron tomography and so forth. 

Among the most prominent compounds intended for thin film productions, stands out the 
Cadmium-Telluride (CdTe) owing to its suitable semiconductor properties such as direct energy 
gap ( E gap = 1.53 eV to 300 K) and high optical absorption coefficient (k 5 x 10 5 cm x ) 

[4] . Applications concerned on CdTe films span the fabrication of solar cells of high efficiency 

[5] as well as of the X-ray, 7 -ray and infrared detectors [ 6 ]. There are extensive studies on 
CdTe growth in branches as diverse as the controlled growth of tetrapod-branched nanocrystals 
[7], self-assembly of quantum dots [8,9], nano-wires [10], microcavities of ultrafast optical 
responses [11], etc. Nevertheless, there are still a few works on CdTe growth dynamics itself 
[12-15], discussing issues as kinetic roughening, growth symmetries, as well as morphological 
aspects like grain size, grain shape and their dependency on deposition temperature 1 , molecular 
flux and thickness. In fact, these are the most important structures and parameters affecting 
elasto-mechanical and electrical film features and, consequently, the efficiency of devices built 

'This term refers to the temperature of the substrate. 
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up from them [1,2, 16, 17]. 

From the theoretical side, the surface (or interface) growth is an interesting out-of- 
equilibrium Statistical Mechanics subject because of scale invariance and Universality emer¬ 
gence, as occurs in thermal fluctuations of equilibrium systems at criticality [18]. The scaling 
invariance implies absence of any characteristic length in the system beyond system size itself. 
In turn, the Universality concept means that systems of different microscopic nature can exhibit 
the same large scale behavior, since they are ruled by interactions sharing dimensionality, sym¬ 
metries and conservation laws [18]. Actually, the Universality idea goes beyond of interface 
studies and has been found in others far-from-equilibrium contexts such as epidemic spreading 
in random networks [19], crackling noise [20] and social dynamics [21]. However, most ef¬ 
forts have been focusing in surface growth due to its ubiquitousness in the nature with examples 
ranging from Physics, Chemistry, Biology to Applied Mathematics and Geology [22], 



Figure 1.1: Fresh snow falling on the window and forming a rough interface. Extracted 
from [22], 


For instance, one can observe the process of snowflakes falling on the window, sliding 
down through it and sticking on the first aggregate they have met. Although the whole system is 
part of our daily experience and seems to be quite “simple”, the generated interface is undoubt¬ 
edly fascinating - see fig. 1.1. One can notice the presence of large voids, branches and a rough 
aspect. 

Complex interfaces (in the sense of formation, growth and dynamic) are found in the 
nature from the bacterial growth, spreading of flame fronts, propagation of fluids in porous 
materials, etc, as well as surfaces made in the lab like those from MBE, eletroconvection and 
other thin-film deposition techniques [22], Figure 1.2 displays a Scanning Electron Microscopy 
(SEM) surface image of Bi film electrochemically grown on a Si substrate. This surface, grown 
on a two-dimensional substrate, is contrasted with that one formed by snowflakes and presents 
a rougher aspect which resembles to fern leaves, an example of macroscopic object exhibiting 
fractal properties. Indeed, interface growth is a close topic to the fractal geometry by Mandel¬ 
brot [23] and a concise link between them is established in the chapter 2. 

In order to describe the dynamic of growing surfaces which are driven by local processes, 
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Figure 1.2: SEM image of Bi thin film grown on Si by electrodeposition technique. 
Courtesy of prof. Rene C. Silva from Universidade Federal de Vigosa, Brazil. 


continuum growth equations can be proposed. These equations does not take into account the 
microscopic nature of systems, but only the underlying symmetries and relevant relaxation pro¬ 
cesses that rule the dynamic at a coarsening-grained level and in a hydrodynamic limit 2 . Under 
appropriate scale transformations, some statistical quantities are kept unchanged and define, this 
way, “critical” exponents. In a successful description, the set of these exponents is related to an 
Universality Class (UC) in which systems sharing dimensionality and large scale behavior can 
be grouped. 

Among a few UCs theoretically predicted, the most remarkable one is that of Kardar- 
Parisi-Zhang (KPZ) [24]. It has became a celebrated UC because: 1- its continuum growth 
equation, which is non-linear, can be mapped in classical equilibrium problems of the Statistical 
Physics [25] and in several mathematical motivated ones [26]. 2- In d = 1 + 1 dimensions 3 , 
important KPZ models as the Single-Step model [27] and the Poly-Nuclear Growth model [28] 
allowed uncovering a parallel between the (rescaled) height distribution of KPZ interfaces and 
the famous Tracy-Widow distributions emerging from the Random Matrix Theory [29,30]. 3 - 
Ten years later, analytical treatments on the basis of the famous Directed Polymer in a Random 
Medium (DPRM) model [25,31] allowed to solve the KPZ(-DPRM) equation [32-37], in the 
meantime that experiments carried out by K. Takeuchi and M. Sano [38-41], regarding the 
growth of turbulent liquid crystals, have confirmed the analytical findings and suggested new 
universal KPZ features. Previous experiments concerning on the low combustion of paper [42— 
45], and a recent experimental realization on the deposit of colloidal particles at edge of water 
drops [46] have also enormously contributed to enoble the KPZ^ =1+1 paradigm. Shortly after 
Takeuchi and Sano experiments [38-41], numerical models belonging to the KPZ class have 
supported and, indeed, gone beyond liquid-crystal results [47-50] filling, finally, the last pieces 
toward a consistent KPZ,; = i + i triumvirate. 

In d = 2 + 1, the KPZ situation, however, is very different from its lower dimensional 

2 Details are found in the chapter 2. 

’This notation means one substrate topological dimension (d s ) plus one growth direction [22]. 
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counterpart. There are not analytical results and almost one knows about the most impor¬ 
tant dimension for technological applications comes from numerical results: the scaling expo¬ 
nents [51,52] and the height [53], squared local roughness [54,55] and maximal relative height 
distributions [56] in the steady-state are older examples. Recently, the universality for height 
distributions in the growth regime has been verified through large-scale simulations [57-60]. 
From the experimental side, evidences of KPZ,/ = 2+1 systems are extremely rare [61,62] and a 
long quest to find out one robust realization confirming the KPZ universality beyond scaling ex¬ 
ponents has persisted until the beginning of 2014 [63,64]. In fact, this poverty of experimental 
evidences has touching all UCs. Several works (see chapter 11 in the Ref. [22] and section 3 in 
the Ref. [65]) have been reporting exponent values that do not match with anyone known UC. 
Inevitably, it has created a deep nuisance in the area, suggesting that the theoretical framework 
as yet is not completed. Advances as the anomalous scaling [66-71] have shown that some sys¬ 
tems are ruled by different exponents at local and global scales and that the analysis of interface 
fluctuations in the Fourier space arc essentials for accessing the true form of the local scaling 
law [70]. However, in these situations, a larger number of exponents must be found in order to 
classify the interface dynamic. It makes the work still harder. 

From an experimental point of view, there are some basic aspects hampering the associ¬ 
ation between a two-dimensional growth and its UC. They are listed in the following: 

i) Difficulty for imaging surfaces and for growing films at long times - Unlike one¬ 
dimensional growing interfaces, it is hard recording the time evolution of two-dimensional 
surfaces during the growth. Thus, in general, ex-situ probe microscopic techniques are 
used for imaging the interface at discrete growth times. Each growth time corresponds to 
an one distinct sample produced in the lab, which, depending on the growth technique and 
on the growth parameters, can take some hours (or even days) to become ready. Because 
of this procedure, two-dimensional growths are, usually, conditioned to the finite-time 
growth, instead to the asymptotic-time growth, where the scaling regime of the true UC is 
expected to emerge. 

ii) Finite-time effects - During the growth, the interface dynamic can suffer a transition, i.e, 
a crossover. Temporal crossovers can be identified measuring the time evolution of some 
correlation function 4 , but even in simulations this is not an easy procedure [66,72-74]. 
For instance, studying the growth of Si02 on Si by CVD, F. Ojeda et al. [61] have found 
the true asymptotic scaling regime only for samples grown at the range of rs_/ 10 2 to 10 3 
min, after two initial crossovers - see fig. 1.3. From this example, it is clear that if the data 
from an experiment corresponds to a crossover region (very difficult to be detected), the 
exponent value extracted from there could not match with any universal value. The same 
can occur when transient anomalous scaling takes place [66]. Both facts are also reasons 
of why so many experimental works have found exponents which are not theoretically 
expected. 

iii) The presence of morphological instabilities - When the surface can be decomposed into 
an array of globules (grains, mounds, etc.), one says there is the presence of morphological 

4 They are thoroughly defined in the section 2.2. 
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time (min) 


Figure 1.3: Example of crossovers ocurring in the two-dimensional growth of Si0 2 
on Si(001) substrates by CVD. The crossovers are identified by different slopes (blue 
to red and red to black) at the W(t) scaling in a log x log plot. The short number 
of samples makes hard to distinguish different regimes. The value 0.26 refers to the 
slope of the black curve and was related to the assymptotic scaling regime. Figure 
extracted and edited from [61]. 


instabilities [75]. This feature induces a characteristic length in the system £, below which 
scaling invariance is awaited to be broken. The presence of £ modifies local scaling of 
correlation functions and has led to many equivocated associations between experimental 
growths and UCs. In particular, non-universal local exponents have been confused with 
critical ones, as very well explained by Oliveira and Aarao Reis [76,77]. These studies 
are addressed in the section 2.3.2. 

iv) Dimensional fragility of the KPZ equation - Concerning on the KPZ class, as theoret¬ 
ically discussed in the ref. [78], a same experimental KPZd=i+i system, when grown in 
d = 2 + 1 dimensions, might has its exponents deviated from those expected for the KPZ 
class. This occurs due to the presence of morphological instabilities and non-locality, 
which are introduced by hand as a perturbation term in the non-linear term of KPZ equa¬ 
tion 5 . M. Nicoli et al. have been calling attention for this dimensional fragility as a further 
obstacle in front of experimental confirmations of KPZ^ =2 +i growth. As extension, one 
can wonder whether this fragility is really a particular feature of the KPZ equation or can 
take place in others non-linear growth equations. 

Guided by all this experience with two-dimensional experimental systems, as well as 
the deep knowledge about the KPZ,/ =1 + ] triumvirate, we have found out the first robust exper¬ 
imental confirmation of the KPZ class in two-dimensional systems, going beyond the standard 
comparison with exponents [63]. Through a very systematic scheme for investigating the UC of 


5 See chapter 3. 





1. Introduction 


6 


growing films, which advances over the comparison of scaling exponents and provide a reliable 
way for circumvent out (partially) the obstacles ii-iv), we have confirmed that the growth of 
CdTe thin films on Si(001) for T = 250 °C belongs to the KPZ class as demonstrated by scaling 
exponents and universal height distributions, universal squared local roughness distributions and 
universal maximal relative height distributions [63]. Moreover, this is the first time that it has 
been experimentally demonstrated the universality of such distributions [63]. 

As the goal of the present work is twofold (theoretically and experimentally-motivated), 
we have also studied the effect of the deposition temperature on this novel two-dimensional 
Kardar-Parisi-Zhang system. We have shown that in the presence of finite-time effects (crossovers 
and anomalous scaling), scaling exponents are not able to point, in a compelling way, the UC 
class of the system, mainly due to inexorable experimental obstacles. However, we have shown 
that this can be circumvented by analysing distributions, which allow us to prove that fluctua¬ 
tions of CdTe interface are described, asymptotically, by the Kardar-Parisi-Zhang equation in a 
broad range of temperature [79]. Finally, one demonstrates that is possible tunning, experimen¬ 
tally, the KPZ non-linearity through the deposition temperature [79]. 

The rest of this dissertation is organized as follows. Chapter 2 makes a link between 
fractals and surface fluctuations, introducing basic statistical tools. Up-to-dated results about 
the specific subjects of each section can be found diluted along the text and, when indicated, in 
the Appendix sections. The chapter 3 is deserved for a short review on the Kardar-Parisi-Zhang 
equation, which includes the mainly ingredients and advances whithin the KPZ theory. In the 
chapter 4, experimental procedures used during this work are carefully described, as well as the 
techniques and growth parameters. Chapter 5 presents the first results of this work and exhibits 
a breaktrough in the context of the two-dimensional KPZ paradigm. On the sequence line, in the 
chapter 6 one distills the effect of the deposition temperature on the CdTe growth dynamic, and 
one explores the richnnes emerging from the CdTe system. Finally, in the chapter 7, conclusions 
and perspectives are drawn with a compreheensive overview of this work, standing out its mainly 
contributions to the context of kinetic roughening surfaces. Appendix sections cover topics 
as the random growth equation, the Villain-Lai-Das-Sarma (VLDS) class, etc., the origin and 
development of anomalous scaling and, finally, an introductory material for young students 
concerning on the basic physics of crystal growth. 



2 


Fractals, Scale Invariance and Universality in 
Interface Growth 


In this chapter, a link between growing surfaces and their descriptions, built on the analogy of 
tools from Statistical-Mechanics of equilibrium phenomena at the criticality is given. Contin¬ 
uum growth equations are discussed and their critical exponents arc extracted. Universality is 
discussed on the theoretical, numerical and experimental approaches. 


2.1 From Fractality to the Family-Vicsek Ansdtze 

A fractal is a complex and irregular 1 object which, under an appropriated scale transformation, 
(any)one of its parts represent it as the whole. Mountains, clouds, tree leaves and coastlines 
are some examples of fractals decorating the nature [23], In a mathematical sense, a fractal is 
said deterministic if a “zoom” in the system always reproduces exactly the whole object. How¬ 
ever, in the nature, one finds statistical fractals. For example, if one compares two snapshots 
from a mountain at different magnifications they do not overlap but, nonetheless, their statistical 
properties are the same [22] . 

A fractal object presents a dilatation symmetry (or homogeneity property) [18]. Statisti¬ 
cal homogeneity is described by the condition: 


f{S ai X 1,?" 2 X2,?" 3 X3, ...) = ?“/(x 1,X 2 ,X 3 , ...), (2.1) 

where the object is formed by a set of points / = /(x), ? is a scale factor and a is called the 
Holder exponent [22]. 

If a\ = c >!2 = = ol (meaning an isotropic scale transformation) satisfies the eq. 2.1, so 

one has self-similarity. Otherwise, for anisotropic transformations, one has self-affinity [22,23]. 

1 In the sense that Euclidean geometry can not describe it. 
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At a critical point of a phase transition, correlation functions 2 behave exactly as eq. 2.1 
[18]. In the Statistical Mechanics language, it implies that the correlation length 3 (0 diverges 
and there is no characteristic length in the system beyond system size itself [18]. For non- 
equilibrium processes such as surface growth, it has been shown that interface fluctuations also 
behave as those at critical point of a phase transition [18,22,65]. In other words, typically, 
interfaces provided by the nature are fractals in the statistical sense. Let h(x, t) be the height of 
an one-dimensional interface at the position x at time t. Assuming homogeneity, as in eq. 2.1, 
then: 


h(x,t) = ? a h(qx,q z f), (2.2) 

where z is the dynamic exponent, and a, in the surface growth context, is named roughness 
exponent [22]. Notice that, the dilation symmetry has also been assumed on the temporal axis, 
which has lead to the z definition. Hence, a and z are characteristic of systems displaying 
spacial and temporal scaling invariance. 

A surface evolving in time is not a deterministic process. The height h{x) at surface is a 
stochastic variable and have a set of possible outcomes {hi, h%, I 13 ...} distributed from a given 
probability P(h). In general, the probability for a given height of the interface take the value h 
is defined as the ratio between the number of occurrences Nf, of this specific event and the total 
number N of heights sampled: 


P(h) = lim Nft/N. (2.3) 

TV—>00 

The probability depends on the range (dh) in which a value between h and h + dh 
is measured. So, is more interesting to define the probability density function (pdf) p(h) = 
dP(h)/dh [18, 80, 81], which must satisfy the normalization requirement: 



dhp(h ) 


= 1. 


(2.4) 


A pdf is built up from its moments or cumulants. The nth -moment ( m n ) of a pdf is: 


/ OO 

dhp(h ) ( h n ), (2.5) 

-OO 

where the brackets mean to take the expectation value of the variable h n . 

The moments are generated through the characteristic function p(k), which is the Fourier 
transform of p(h) [80]. From the logarithm of p(/,•), one defines the nth -cumulant ((h n ) c ): 


lnp(fc) = f;LiO0"}c. (2.6) 

■' to ! 

n=l 

2 A measurement of how local fluctuations in one part of the system affects those in the other ones. They are 
defined explicitly in the next section. 

’The length over the which one local part of the system affects the other ones. 
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Figure 2.1: Comparison between different distributions having the same expectation 
value, which is indicated by the solid red line. A curve presenting S < 0 can be obtained 
performing a parity transformation on the green curve. 


By expanding the logarithm in the p(k) definition and comparing therm by therm with 
the expansion of p (k) in powers of k, one reaches to the relation between the moments and the 
cumulants [80]. For the first four cumulants, one has, eq. 2.7: 


(h)c = (h >, 

(h\ = (h 2 ) - (h) 2 , 

(h 3 ) c = (h 3 ) — 2>(h 2 )(h) + 2 (h) 3 , 

(h 4 ) c = (h 4 ) - 4 (h 3 ){h) - 3 (h 2 ) 2 + 12 (h 2 ){h) 2 - 6(h) 4 . 

The first cumulant is called mean and the second one is variance. 


(2.7) 


Other important quantities are dimensionless cumulant ratios. In particular, the Skewness 
S (eq. 2.8) is an indicative of pdf asymmetry [80]. If the expectation value is larger than the 
more probable one, then S > 0; otherwise, S < 0. The Kurtosis K (eq. 2.8) provide information 
about the “weight” of pdf tails. A symmetric pdf with K > 0 presents a peak sharper than that 
from the Gaussian and its tails take longer to fall down. The opposite occurs for K < 0. By 
definition, the Gaussian pdf has S = K = 0. A representation comparing distributions with the 
same expectation value for different S and K values is shown in the fig. 2.1. 


S = (h 3 ) c /[(h 2 ) c } 3/2 and K = (h 4 ) c /[(h 2 ) c } 2 . (2.8) 

The global squared roughness ( w 2 ) of an interface is defined as the variance of heights 
composing it 4 . Interestingly, the roughness (or width) is a very important variable from both 

4 In surfaces with translational symmetry, this is also the variance of p(h). 
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experimental and theoretical point of view. The electrical conductivity of some thin films, for 
example, is strongly reduced as rougher is their surface [1-3]. The importance from the theoret¬ 
ical side will be clear soon. Consider an one-dimensional interface of size L. The roughness of 
such interface at time t reads: 

w(L, t ) = [(h(x, t ) 2 ) - (h(x, f)) 2 ] i/2 , (2.9) 

where the (...) refers to a spatial average over the whole system of size L. 

Now, we can return to the hypothesis made in eq. 2.2. Inserting this equation in the eq. 
2.9 and performing algebraic manipulations 5 one finds: 

w(L,t) =t a ^ z f{L/t^ z ), (2.10) 

this equation precede a canonical ansdtze in the kinetic roughening theory. Defining the growth 
exponent (5 = a/z, the Family-Vicsek (FV) dynamic scaling ansdtze predicts that / asymptot¬ 
ically behaves as f(v) ~ u a for u <C 1 and ~ const for u A> 1 [82]. Leading the ansdtze for 
eq. 2.10, finally, one obtains: 


w(L, t) ~ 


tP, for 

L Q , for f 1 / 2 » L. 


( 2 . 11 ) 


Roughness behavior is sketched in the fig. 2.2. The system roughens as while the 
correlations are spreading through it. This is called growth regime. At the time t x , the correla¬ 
tion length becomes of the same order of L and a distinct regime is reached: roughness stops 
growing in time and turns to depend solely on L. This is the steady state or saturation regime. 
Comparing what eq. 2.11 tell us, it is clear to associate £|| ~ f 1 / 2 , where £|| is the parallel 
correlation length. In inset of fig. 2.2, the FV ansdtze is tested, showing remarkable collapse for 
the curves. We notice that the roughness presents a power-law dependence in space and time, 
as occurs with the correlation functions in equilibrium critical phenomena. The parallel goes 
further and, indeed, the “critical” exponents (a and z) do not depend on microscopic details of 
the system under investigation, i.e, there is universality in fluctuations of growing interfaces. 

Once the saturation regime depends on t 1 / 2 be of the same order of L, in experimental 
situations where L is much larger than the characteristic size of particles constituting the inter¬ 
face, the time required to the system gets into the steady state is hardly achieved. As far as we 
know, only the evolution of profiles (d = 1 + 1) yielded by slow combustion of paper [42] and 
the growth (d = 2 + 1) of Si02 films by CVD (after 2 days of deposition) have experimentally 
achieved the stationary state [61]. 

The FV ansdtze has been confirmed in several examples of surface growth such as the 
propagation of fluid flow in a porous medium [22], paper wetting [22], the growth of turbulent 
liquid crystals [38-40], the slow combustion of paper sheets [42,43] and so on [83]. However, 
the FV scaling is not the most general one. Indeed, FV ansdtze fails for describing local scaling 
of growth process with anomalous roughening. A detailed description of anomalous scaling is 
let to the appendix section B. 

5 One must assume q z t = 1. This can be done because <; is just an arbitrary factor. We can change its value until 
obtain the result q z t = 1. 
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Figure 2.2: Typical time roughness plot (log x log) during an interface growth for 
different substrate sizes. The curves refer to substrate of size: black circles (L), orange 
left triangles (2L), blue up triangles (4L), green diamonds (8L) and red squares (16L). 
t x indicates the characteristic time in which £|| « substrate size. Inset confirms the FV 
ansatze, described in eq. 2.11 using the appropriate values for f3 and 2 . 


2.2 Correlation Functions 

Correlation functions (CF) play a central role in systems exhibiting criticality (out-of or at equi¬ 
librium) because they provide a measurement of the correlation length. They are formed by an 
operation (sum, product, etc.) involving suitable quantities describing the microscopic state of 
the system 6 , which are separated by a distance of l. Taking an interface of linear size L for in¬ 
stance, one can describe it by its height field {K}l or even by its slope field {V/i} l (of course, in 
an appropriate coarsening grained level). A CF for an interface growing can, hence, be written 
as: 


C h (\l\,t) = ([ h(x,t ) - h(x + l,t)} 2 ). (2.12) 

If the system presents translational symmetry, any CF will depend solely on the magni¬ 
tude of l [18], what is not true for anisotropic systems (see [84] and ref. therein). Furthermore, 
as discussed at the beginning of this chapter, a CF displays a dilation symmetry at the criticality. 
Hence, inserting eq. 2.2 in eq. 2.12 and performing algebraic manipulations one reaches to: 

C h (l,t) = t 2 a ' z f(l/t 1 / z ), 

6 Examples are the local magnetization for a spin lattice, the local height for an interface, etc. 


(2.13) 
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where f(u) is the scaling function that, for while, obeys the FV ansatze. Ch is often called 
height-difference correlation function. 

Notice that this is the same scaling form obtained for the roughness in the eq. 2.10 with 
/3 replaced by 2/3 and L by l. So, at first note, the example showed in the fig. 2.2 is shared by 
Ch(t), but with the growth regime evolving as t 2 ' 1 . In the same way, one can replace L by I in the 
eq. 2.9 to obtain a scaling form for the local roughness [vji oc (1, t)], which is equally a correlation 
function 7 . From an experimental point of view, it is unpractical changing L of a system in order 
to extract a. Rather, one often uses local measurements, spanning boxes of lateral length l in 
the interval [0, L\ and obtaining the exponent from the hypothesis CJ ( l , t) Wlocijff) 

Fig. 2.3 shows a typical behavior of wi oc for interfaces in the growth regime. 



/ 

Figure 2.3: Typical local roughness plot (log x log) for an interface grown at different 
times. The curves refer to growth times: black circles (t), orange left triangles (2t), 
blue up triangles (4t), green diamonds (8t) and red squares (16t). Vertical dashed 
lines indicate the £|| measured. Inset confirms the FV ansatze, described in eq. 2.11 
with L replaced by l, and using the appropriate values for a and 0 . 


Other important CF is the spatial covariance of heights (C s ), defined as: 

C s (|l|,f) = ([h(x,t)h(x + l,t)\) - (h) 2 , (2.14) 

where it is straightforward to show that Ch = 2 w 2 oc — 2 C s . 

The covariance C s (l,t ) has been computed, in particular, for one-dimensional models 
belonging to the Kardar-Parisi-Zhang (KPZ) Universality Class (UC) because they are universal 
and given by the covariance of Airy processes [28, 38, 59] (see fig. 5 in the ref. [38]). Very 


7 In the absence of anomalous roughening, the w; oc scaling follows the FV ansatze as in the eq. 2.10. 













2. Fractals, Scale Invariance and Universality in Non-Equilibrium Processes 


13 


recently, T. Halpin-Healy and G. Palazantzas have found numerically and confirmed experimen¬ 
tally the universality of the rescaled C s (kpz) also in 2 + 1 dimensions [64]. 

Likewise, one can define slope-slope correlation functions using Vh, instead of h, in the 
eq. 2.12 and eq. 2.14. Indeed, several experimental studies [15,63,85] have used the slope-slope 
covariance (eq. 2.15) in order to obtain an estimative of £||. 

r(|l|,i) = ([V/i(x, f)V/r(x + l, t)]). (2.15) 

An estimative of £| | can be done measuring either the first zero or the first minimum of 
the curve [17]. Figure 2.4 shows a plot of a typical behavior of the slope-slope covariance as 
function of the distance 1 and time t. The procedure used for estimating £|| is also indicated, as 
well as the plot of £|| as function of the growth time, from where the exponent 1 /z can be found. 



/ 

Figure 2.4: Typical normalized slope-slope covariance as function of the distance I 
and the time t in a linear x linear plot. Curve colors refer to growth times: black (t), 
orange (2t), blue (4t), green (8t) and red (16t). Magent dashed line sets r(|l|,t)/T(0,t) 
equal to 0. Vertical dashed lines indicates the position of the estimated correlation 
length. The inset shows £||, as function of time in a log x log plot, from where the 
exponent l/z can be extracted. 
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2.3 Continuum Equations and Universality Classes 

Let an interface be described by its height field {h(x,io)}, in a appropriate coarsening grained 
level, and consider there is a known driving force, i.e., there is a far-from-equilibrium situation. 
An important question concerning on it could be: how does one can describe the height-field 
evolution of the interface? A priori, the answer does not appear to be reachable because it seems 
that there is a lack of informations such as: i) the kind of interface (biological, physical, chemical 
reaction front, etc.) which is being dealt and ii) the specific interactions ruling the dynamics at 
the microscopic level. 

Indeed, an approach to the question from this point of view makes the problem practi¬ 
cally intractable. Nevertheless, looking at the macroscopic scales of evolving interfaces, one can 
glimpse that the collective behavior of them shares many underlying similarities, which does 
not depend, in fact, of i) and ii). Hence, based on this one could perform a long-wavelength 
description of the system. Advancing in these thoughts, one can argue that, beyond of the 
long-wavelength hypothesis, the long-time regime is also a necessary condition, since transient 
(finite-time) behaviors should be avoided. This particular limit of long wavelengths and long 
times is called hydrodynamic limit by analogy to the Navier-Stokes equation for a fluid of parti¬ 
cles [18]. 

Guided by these ideas, one can build general equations for describing the height-field 
evolution of growing surfaces, with the form: 

dth(x, t) = F + 0(x, {h},t) + rj, (2.16) 

with F being the average number of particles per unit time arriving at the interface 8 , the called 
driving force. Now, one must account that this arriving process is stochastic, and the therm rj is 
inserted in the equation for capturing this feature. Regarding 77 , one has: 

• If the active zone (interface front) advances onto an inhomogeneous medium as a porous 
substrate or a paper sheet, the relevant noise in the process is static and changes point 
to point in the medium. This is called of quenched noise, once rj = p(x, h ). If only 
the quenched noise is present, the growth of the interface has a deterministic evolution. 
However, the thermal noise, which is always present in experiments, destroys this deter¬ 
minism. 

• Beyond of thermal noise, when an interface grows by receiving particles from an external 
flux, the so called shot noise is also present and plays the crucial role on the dynamic. This 
noise comes from the inherent randomness occurring in the deposition process. Assuming 
there is no preferential area onto the substrate for receiving molecular flux, one has that 
(t/(x, t )) = 0 with the spatial-temporal covariance given by: 

(r/(x, t)rj(x , ,t')) = 2 D5 da (x — x’)5(t — t'), (2.17) 

where \fD is the amplitude of the white noise. 

8 Assume the particle size being equal the unit, in order to satisfy dimensional analysis. 
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Now we turn the attention to the forms which the functional 0(x, {h},t) can assume sup¬ 
posing that the dynamic of the interface is ruled by local processes. As an underlying hypothesis, 
descriptions of a physic fact can not depend on the origin of its observation a = (to, x, ho). 
Thus, spatial and temporal translations must be satisfied by the eq. 2.16. It rules out from the 
functional explicit terms involving t n , x m or h m , where m G M*. Moreover, as the growth does 
not make distinction between “right”- and “left-handed”, the eq. 2.16 must also be invariant 
under spatial parity transformations with respect to the x axis. These hypothesis reduce the 
allowed terms in 0 to combinations of even derivatives such as (V 2n h)('Vh) 2p , with n,p E N. 
In a general picture, our considerations until here have led us to consider: 

d t h(x, t) = F + aiV 2 /i + a 2 V 4 /i + ... + h (V/i ) 2 + 6 2 (V/i ) 4 + ... 

+c n (V 2 /r)(V/i ) 2 + ... + c np (V 2n h)(\/h) 2p + rj, (2.18) 

where a*, b, and c np are appropriate quantities making the eq. dimensionally consistent. 

As we are interested in the hydrodynamic limit, derivatives of higher order are irrelevant 
to the asymptotic scaling behavior (see ref. [22] pag. 49). Thus, the simplest general equation 
involving these terms reads: 

dth(x, t) = F + ai V 2 h + h (V/r ) 2 + cu (V 2 /i) (V/t ) 2 + 77 , (2.19) 

As we shall see, important continuum growth equations arc encoded in the equation 
2.19. Nevertheless, there arc also others important growth equations that call for a physically- 
motivated term which makes appear higher derivatives of h. In the following section we discuss 
this subject in details. 

2.3.1 Edwards-Wilkinson and the Linear-MBE equation 

The Edwards-Wilkinson (EW) equation (eq. 2.20) was proposed in 1982 for describing the 
sedimentation of granular particles [ 86 ]. The equation preserves parity symmetry in the growth 
direction with respect to the mean height, leading the height pdf to has S = 0. This particular 
symmetry rules out even powers of h such as (V/i) 2n and, in accordance with the general eq. 
2.19 9 , it reads: 


dth(x, t ) = uV 2 h + r](x, t ), (2.20) 

where v = a\ is a quantity of dimension [m 2 /s] in the S.I. convention. 

The growth average velocity of the interface v = (dth) = Ft, once (V n /i) vanishes 
for periodic bound conditions. The eq. above is written at the referential of the mean height, 
avoiding the explicit dependence on the flux term and setting v = 0. 

Physically, the laplacian term acts as a conservative smoothing mechanism redistributing 
the irregularities on the interface, while maintaining the average height unchanged (see pag. 50 
in the ref. [22] for a geometric interpretation). 

9 One can show that the term (V 2 /i)(V/i) 2 under renormalization is irrelevant compared with (V 2 /i). See Ref. 
[22] pag. 49. 
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Due to the linear character of this equation, critical exponents can be found straight¬ 
forwardly by rescaling or by Fourier transform methods. As a prototype case, we shall find 
exponent values by applying rescaling tools on a general linear equation obtained replacing 2 
by 2 n on the nabla operator in the eq. 2.20. The method by Fourier transform can be found in 
details in the pag. 173 of the Ref. [83]. 

Rescaling: Suppose a rescaling as x —> gx, t —> <, z t and h —t g a h. Inserting in the eq. 
2.20 with 2 replaced by 2 n on the nabla, one has: 


= vq a - 2n V 2n h + r/( ? x, g z t), (2.21) 

where, using the noise covariance definition (eq. 2.17) and the delta function properties 5 ds (<,x) = 
S~ ds 5(x), one can rewrite the rescaled noise in the expression above by ^~ z ~ d3 ^ 2 r](x, t). Now, 
assuming scale invariance, one finds the critical exponents (valid for n > 1): 


a = 


2 n — d s 
2 


and z = 2 n. 


( 2 . 22 ) 


The EW eq. is recovered for n = 1, which yields a = (2— d s )/2 and z = 2. In particular, 
for d s = 2, one obtains a = 0 and (5 = 0 meaning that in the eq. 2.11 the roughness exhibit a 
logarithm dependence on t and on L. The set of exponents a = a(d s ) and z = 2 compose the so 
called Edwards-Wilkinson Universality Class, well described in the hydrodynamic limit by the 
EW equation. A famous numerical model belonging to the EW class is the Random Deposition 
with Surface Relaxation (RDSR) proposed in 1986 by Family [87], where the deposited particles 
are allowed to diffuse at surface until reaching the local lowest height. Experimental evidences 
of surfaces belonging to the EW class, in turn, are very rare. As far as we are concerned, 
the EW universality has only been found in the growth of W multilayers on Si by magnetron 
sputtering [88,89] and in the sedimentation process of Si02 nanospheres [90]. 

Returning to eq. 2.22, if one sets n = 2 one obtains the exponents for the famous growth 
equation known as linear-MBE equation (eq. 2.23), firstly proposed by Wolf and Villain [72] 
and Das Sarma and Tamborenea [91] for describing the growth of surfaces in which diffusion is 
the relevant growth mechanism. 


d t h(x, t ) = -K d V 4 h + rj(x, t), (2.23) 

K,i accounts for the strength of the diffusion and has dimension of [m 4 /s]. 

A derivation of eq. 2.23 from a conservation law is let to the appendix section A.2. 
Indeed, the deterministic form of eq. 2.23 was known (and was solved) since a long time ago 
by Herring [92] and Mullins [93] considering the effect of the scale on sintering phenomena 
and the development of thermal grooves, respectively. Due this, often the stochastic growth 
equation is called Mullins-Herring equation, instead linear-MBE. However, in this dissertation 
we will refer to the eq. 2.23 as linear-MBE equation, while its related universality class is being 
called Mullins-Herring (MH) class. 

The critical exponents constituting the MH class are (from eq. 2.22, with n = 2): 

4 - d s 


a = 


2 


and z = 4. 


(2.24) 
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Theses exponents were firstly calculated by Wolf and Villain as intention of describing 
the WV model [72], although, shortly after was found that, actually, the model does not belong 
to the MH class [74]. Numerical models capturing the diffusion mechanism and truly belonging 
to the MH class carries the name of Das Sarma, Tamborenea, Ghaisas and J. Kim [94,95]. A 
classical review paper on the models can be found in the ref. [96], while a discussion on the WV 
model is described in the chapter 15 of the ref. [22]. 

On the experimental side, evidences of MH class have been found in the growth of Si 
on Si(lll) by MBE for T = 275 °C and deposition rate of 7 bilayers/min [97] [scaling was 
obtained through STM images]; thermal evaporation of amorphous Si on Si(lll) substrates 
at deposition rate of 0.8 ± 0.2 A/s [98] [AFM]; the sputter-deposition growth of Pt on glass 
at deposition rate of 6 A/s and with the normal substrate aligned about 45° with the target 
surface normal [99] [STM]; in fluctuations of intra-grain domains of gold electrodeposits grown 
at 100 nm/s [100] [STM] and in inter-grain fluctuations of LiCoO x thin films grown by rf 
sputtering, after annealing process [101] [AFM], as well as in the circular growth of cultivated 
brain tumor [102] [optical microscope]. 

2.3.2 The Linear Deposition-Desorption-Diffusion equation and Grain Models from 
Oliveira and Reis 

At this point, we could insert together the diffusion and the desorption mechanisms in a contin¬ 
uum growth equation. This makes sense, once specific MBE conditions (very high deposition 
temperature and/or low supersaturation) can lead the growth to depend sensitively on both pro¬ 
cesses. Deposition occurs solely if there is a supersaturation ruling the vapor-to-solid phase 
transition, it means that the difference of chemical potential between the vapor (//„) and above 
the solid phase (/r(x, /)) is positive. However, when the supersaturation is low [ (p v — /;) > 0], 
thermal fluctuations make possible, during a characteristic time, a particle to escape from the 
interface and to move toward the vapor phase. In these conditions, and including diffusion, one 
has: 


d t h(x, t) = F - B[(/j, v - (j,)] - K d \7 A h + rj(x, t ), (2.25) 

with B having dimension of inverse of linear momentum. 

Inserting fi(x, t) oc — V 2 /i, one finds the general deposition-desorption-diffusion equa¬ 
tion: 


dth(x, t) = F + vS7 2 h — K d V 4 h + B(i v + 77 (x, t), (2.26) 

where the ratio {K ( ]/u)^' 2 has dimension of length and, thus, define a characteristic length in 
the system, £. 

A trivial rescaling in the eq. 2.26 provides: 

q a ~ z d t h = i/c°- 2 V 2 h - Kd^-^h + BfjL v + ? - z - d v ? (x, t), (2.27) 

where for long-wavelength fluctuations (c — y 00, l ^$> Q, the laplacian dominates the growth 
and the exponents are consistent with the EW class. At short-length scales (c -» 0, l <C Q, 
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however, diffusion mechanism overcomes the laplacian effect and the growth is dictated by the 
linear-MBE equation. 

This behavior has been confirmed, for instance, in copper electrodeposition in the pres¬ 
ence of l,3-diethyl-2-thiorea, an organic additive with concentration x (0.3 < x < 0.4 mM), 
for low current density j = 0.02 A cm~ 2 [103]. In this work, the authors have found a and 
/3 exponents for both regimes, strongly confirming the MH-EW crossover. Notwithstanding, in 
many studies on the growth of thin films the local roughness is the only curve analyzed, from 
where a is usually extracted. In particular, when a « 1.00 a linear-diffusion-dominated dy¬ 
namic is often suggested, but Oliveira and Reis have demonstrated that this procedure is fail in 
most of cases [76,77]. 

Studying the effect of grains at surface on the local roughness scale, Oliveira and Reis 
have shown that the general uy oc (r, t) 10 plot exhibits two crossovers - see and follow the fig. 
2.5. The first regime happens for r <C r c , where r c is the average grain size and is dictated by 
the exponent a±. Length scales between grains r r <C r <C £ define the second regime growing 
as wi oc ~ r" 2 . The second crossover separates correlated from non-correlated regions of the 
grainy surface. 



Figure 2.5: Grain of (a) semi-elliptical, (b) conical and (c) pyramidal shapes, (d) 
General local roughness curve (log x log) for grainy surfaces. Fig. extracted and 
edited from Oliveira and Reis [76,77], 


Interestingly, it has been demonstrated that a\ has a dependence on the grain geometric 
shape, on the local correlation function and on the procedure to calculate root-mean-square aver¬ 
ages. It implies that a\ can not be related to a critical exponent in the sense of capturing univer¬ 
sal fluctuations at surface. The geometric interpretation is corroborated with the local roughness 
curve, where the cvi value can change from 0.75 to 1.00 for surfaces composed by pyramidal and 
flat top grains, respectively. The 012 exponent, however, keeps constant at the expected UC value 
related to the model. We recommend the reader take a look at the original papers from Olivera 

in We are preserving the original notation of the paper [77], 
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and Reis [76, 77]. Comparison with experiments has been carried out and, just to few some 
experimental results matching with Olivera and Reis predictions, one finds: the spray pyrolysis 
growth of ZnO films [104], which gives 0.94 <C «i <C 0.97 for high flow rates; the electrode¬ 
position of cooper by Mendez et al. [105], presenting a\ = 0.87 ± 0.06; the sputtering of niquel 
oxided films (cei = 0.70) [106]; the growth of bilayers of poly(allylamine hydrochloride) and 
a side-chain-substituted azobenzene copolymer (Ma-co-DR13), after deposition of 10 or 20 bi¬ 
layers, giving a\ ~ 0.80 [107]; and the surfaces of Langmuir-Blodgett films of polyaniline and 
neutral biphosphinic ruthenium complex, which yields 0.66 Car C 0.82 [108]. 




The Kardar-Parisi-Zhang Universality Class: A Brief 
Historical and State-of-the-Art 


This chapter is a natural sequence of the previous one, but due to the richness of the Kardar- 
Parisi-Zhang equation an entire chapter is deserved to the discussions that one follows. A basic 
narrative on the motivation, development and current status of this paradigmatic continuum 
growth equation is presented. We do not have the intention to englobe in this text all what 
the equation has been demonstrating to offer. Instead, we focus on its scaling properties, on 
important mappings and on its universal distributions. 


3.1 Scaling, Mappings, Height Distributions 

The Kardar-Parisi-Zhang (KPZ) equation was proposed in 1986 for describing growing inter¬ 
faces in which growth in the local normal direction plays a decisive role in the asymptotic 
dynamics [24]. Based on universality ideas, as well as on existing growth models as the Eden 
model for cell colony formation [109], and the ballistic deposition model for colloidal aggre¬ 
gates [82, 110], KPZ proposed the simplest nonlinear stochastic equation for describing the 
height-field dynamic of such interfaces: 

^ = u\7 2 h + ^(Vfe) 2 + rj(x,t), (3.1) 

v represents a “surface tension” due to the geometric interpretation of laplacian term (see ref. 
[22]), whereas the nonlinear term accounts for the growth in the local normal direction and rj is 
the noise. 

As indicated in the inset of fig. 3.1, when an interface grows laterally, the increment 
along the h axis (dh) and the local normal velocity v are related by the Pythagoras’ theorem: 


20 





3. Kardar-Parisi-Zhang Universality Class 


21 




Figure 3.1: Successive profiles for a deterministic KPZ growth with a rough initial 
condition. Inset indicates how local normal growth occurs. Extracted and edited from 
[24]. 


5h = [(v5t) 2 + (vdtVh) 2 ] 1 / 2 . Inserting the small slope approximation (V/i <C l) 1 , one can 
perform an expansion at the right side of the equation leading to: ^ ~ v + (v/2)(Wh) 2 + ... 
. Eq. 3.1 is obtained after a transformation to the moving frame (h —> h + vt) and after the 
relaxation mechanism is inserted. Nonlinear higher-order terms are disconsidered because they 
go to zero faster than (V/i) 2 in the hydrodynamic limit, as already discussed in chapter 2. 

Unlike linear growth equations, the mean interface velocity of the KPZ equation takes 
the form: v = Ft + | d da x((S7h) 2 ), which is different from zero even when F = 0. A 

accounts for this “excess of velocity”. In particular, if A = 0 one obtains the EW equation (eq. 
2.20) and if, additionally v = 0, one finds the Random Growth equation (see appendix section 
A.l). 

The KPZ equation can be easily solved in its deterministic form (r/(x. t) = 0), which 
leads to interfaces composed by paraboloid segments resembling dendrites [24], see fig. 3.1 for 
a typical one-dimensional growth pattern. On the other hand, the stochastic equation (eq. 3.1) 
has been resisting to a complete analytical handling. Great part of the advances on the solution 
of the KPZ equation has moved on owing to the arsenal of powerful tools that mathematicians 
has brought to the field. For instance, analytical models that belong to the KPZ equation as 
the Single-Step model [27] and the Poly-Nuclear-Growth (PNG) model [28, 111] have, togheter 
with, the Totally Asymmetric Exclusion Process (TASEP) [112] and the Directed Polymer in a 
Random Potential (DPRM) [31] (which can mapped onto the height field of a KPZ interface), 
played a crucial role within the KPZ theory. For this last, at the beginning of this decade, analyt¬ 
ical treatments have been achieved for d = 1 + 1 dimensions for different geometries, which in 
the surface growth context refers to a curved [32-35], flat [36] and stationary [37] growth Initial 
Conditions (IC). Solutions of higher dimensions, however, are in a fog of analytical hopes. 

Regarding to the KPZ scaling, if one applies, naively, the trivial rescaling used in the 

1 This approximation is not required in others derivations, whose large scale limit take form of the KPZ equation 

[ 22 ], 
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section 2.3.1, one obtains three self-inconsistent scaling relations. Arguing that the nonlinear 
term should dominate the growth in the hydrodynamic limit and applying the rescaling again 
(without the laplacian term) one finds a = (2 — d s )/ 3 and /3 = (2 — d s )/{ 4 + d s ), whose 
predictions are quite different from numerical results [25,47,48,51, 113, 114]. Indeed, this 
procedure is wrong because terms like u, A and D do not renormalize independently - being 
coupled to each other [22]. The correct prediction for KPZ exponents can be achieved using 
Renormalization-Group techniques and mappings to other problems. Through the u = —d x h 
transformation, the KPZ equation with A = 1 is mapped in the Burgers’ equation with noise 
[115], describing the vorticity-free velocity field of a stirred fluid. Additionally, the fact that 
A = 1 in the Burgers’ equation points out that we must preserve the scale invariance on the 
nonlinear term in the KPZ context. This reasoning provides the hyper-scaling relation: 

a + z = 2. (3.2) 

This relation is a consequence of the Burguer equation to be invariant under a galilean trans¬ 
formation, u(x, t) —> uq + u(x — vot, t), which in turn, leads the KPZ equation to be invariant 
under a tilting transformation by an angle e: 

h! —> h + ex; x' —> x — eXt; t' —> t. (3.3) 

The fluctuation-dissipation theorem [22] reveals that, at the steady state (for d = 1 + 1), 
V h follows a Gaussian distribution likewise the position of a particle in a Brownian motion and 
brings out a = 1/2. Together with eq. 3.2, this implies z = 3/2 and f3 = 1/3, in agreement 
with previous [25, 113] and recent [47,48,50, 116] numerical models belonging to the KPZ 
universality class. 

The concept of “universality beyond exponents” has initiated in the classical paper from 
Krug, Meakin and Halpin-Healy, back to 1992’s [113]. In that work, by using the so-called 
“Krug-Meakin toolbox” [117, 118], it was possible to write universal amplitudes in terms of 
model parameters, namely, A = f(I)/i/) 2 and A, which can be easily obtained in simulations 
[48,50,64,113]. Almost 10 years after this imporant result, Johansson studied a model for which 
several probabilistic interpretations can be given [27]. Among them, an interpretation in terms 
of the one-dimensional TASEP (in turn mapped onto the interface fluctuations of the Single-Step 
model) allowed to show that the pdf of a random amplitude (\) related to the height held of the 
Single-Step model (eq. 3.4) is the celebrated Tracy-Widom (TW) distribution [29,30], emerging 
from the Random Matrix Theory context. 

h{t) ~ Voot + sign(X)(Tt) l3 x- ( 3 - 4 ) 

Here is the asymptotic velocity of the interface [ = limt^^oo (d/h) ], sign(X) is the 

signal function, T = aA 1,/Q: |A|, with a being a constant and A = f(D/u) ind= 1 + 1 dimensions 
(see Ref. [113]). 

Shortly after, Prahofer and Spohn [28, 111] unveiled a dependence of p(y) on the geom¬ 
etry (IC) of the growth. By using the mapping of the PNG model onto random permutations, 
they showed that h(x, t ) corresponds to the length of the longest increasing subsequence of such 

2 This relation is valid only for d = 1 + 1 dimensions. In higher dimensions A can also be function of A [113]. 
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X 


Figure 3.2: Universal distributions obtained by Prahofer and Spohn in the Ref. [28]. 
From left to rigth: universal pdfs (solid lines) for x under curved, flat and stationary 
self-similar growth, respectively. Symbols refers to simulations of the PNG model. 
Figure extracted from the ref. [28]. 


permutation, which in turn, is distributed according to: 1) the Gaussian Orthogonal Ensemble 
(GOE) * 3 TW distribution for a flat IC; 2) the Gaussian Unitary Ensemble (GUE) 4 TW, if the 
growth starts from a seed and develops a curved interface; and 3) the Baik-Rains Fq limiting 
distribution [119] for a steady-state IC. Figure 3.2 shows p(y) for different growth geometries 
and for stationary initial configuration calculated from the PNG model and compared to the 
respectives GUE, GOE and Fq distributions [28]. 

A decade later the results of Prahofer and Spohn [28, 111], analytical solutions on the 
one-dimensional KPZ(-DPRM) equation, already cited above, have confirmed the limiting dis¬ 
tributions of x as GUE for curved growth [32-35], GOE for flat growth [36] and Fq for the 
stationary initial condition [37]. Moreover, finite-time corrections were found, suggesting new 
universal KPZ features such as a shift in the mean converging to GOE (GUE) values as t~ r ‘. It 
was also demonstrated that the limiting processes, giving the one-dimensional heigth profiles of 
a KPZ interface, are dictated by the Airyi [120-122] and Airy 2 [122,123] processes for flat, and 
curved growing, respectively. 

Simultaneously, experiments of unprecedented statistics on turbulent liquid crystals per¬ 
formed by K. Takeuchi and M. Sano [38-41] have confirmed carefully great part of the predic¬ 
tions cited above and have given to the KPZd = i + i theory a reliable reality beyond mathematical, 
numerical, and few experimental realizations constrained to exponent results [124]. By using 
the Krug-Meaking toolbox for unearthing non-universal parameters A and A, K. Takeuchi and 
M. Sano calculated the pdf of the \ variable, defined according to the KPZ ansatze (eq. 3.4) as 

’This means, the distribution of the largest eigenvalues of orthogonal matrice ensembles, whose elements are 

distributed from a Gaussian. 

4 Same as before, but now the matrices are unitary. 
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Figure 3.3: Results for turbulent liquid crystals extracted and edited from Takeuchi 
[39]. (a) Histogram of rescaled local height x = h - v aa t/(Tt) 1/3 . The blue and red 
solid symbols show the histograms for the circular and flat interfaces, respectively. 
Dashed and dotted curves show GUE and GOE TW distributions, (b) How x cumu- 
lants approach to the GUE values. Inset shows the mean shift decaying as t~ 1/3 . The 
same behavior is shared by the flat case (not shown here). 


X = (h — v OQ t)/( Ti) 1 / 3 . In the fig. 3.3(a) one can see a comparison between the experimen¬ 
tal p(x) obtained in that work for both curved and flat cases and the GUE and GOE distribu¬ 
tions [39]. Wonderfully, a great accordance was obtained apart from a slight shift at mean of the 
distributions. In fact, in the figure 3.3(b) the shift at the mean is also calculated, suggesting that 
it decays as ~ U l /3 , while the shift for higher cunrulants vanishes quickly. 

Remarkable experiments on colloidal particles deposited at the edge of evaporating drops 
were also able to confirm the “KPZ universality beyond exponents”, where curved deposits of 
particles slightly anisotropic following the GUE-TW distribution were found [46]. Analytical 
and experimental results have been reinforced by numerical ones [47-50], which have filled the 
last pieces toward a robust and consistent KPZ,/ =1 + 1 triumvirate. 

The KPZd= 2 +i situation is very contrasting with its one-dimensional counterpart and 
almost ah one knows about the most important dimension for applications has come from sim¬ 
ulations and, very recently, from some remarkable experimental efforts [63,64]. Best estimates 
for scaling exponents indicate that [51,52,57]: 

a « 0.393 (5 k, 0.242 1/* » 1.607. (3.5) 

The universality of dimensionless cumulant ratios of h, namely, the skewness and kur- 
tosis (eq. 2.8) were firstly calculated numerically at the stationary state, where their values 
were proved to be universal (see Ref. [53] and references therein). The universality in the 
growth regime, although glimpsed in the Ref. [133], only was convincingly demonstrated in 
the year of 2012. Through large-scale simulations, T. Halpin-Healy [57,58] and T. J. Oliveira et 
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al. [59] have uncovered the existence of geometry-dependent KPZ universal height distributions 
at the growth regime, higher dimensional GOE- and GUE-TW counterparts, lying in the heart of 
KPZrf = 2 +i universality. Although the exact forms of these distributions are not known, Oliveira 
et al. have demonstrated that rescaled height pdf’s can be well fitted by generalized Gumbel 
distributions [125] (see the definition in the section 3.3) with parameters m = 6 and m = 9.5 for 
flat and curved cases, respectively. Moreover, as in d = 1 + 1 dimensions, results from the T. 
J. Oliveira et al. study [59] have supported a generalization of the KPZd= 2 +i ansatze (eq. 3.4) 
inserting appropriate finite-time corrections. The final KPZ ansatze reads: 

h(t) = Voo t + sigxiTtfx + V P + C P t~ lp + (3-6) 

where r) p , f p and 7 p are non-universal parameters [60, 126]. The values for these model- 
dependent parameters can be found in the references [47,48,50,60, 113] for d = 1 + 1, and 
[57-60] for d = 2 + 1. Deserving much more attention, the universal-KPZ values for the cumu- 
lants of height distributions are grouped in the tables 3.1 and 3.2. 



GOE 

GUE 

Baik-Rains 

(x)c 

-0.76007 

-1.77109 

0 

(X 2 )c 

0.63805 

0.81320 

1.15039 

ISI 

0.2935 

0.2241 

0.35941 

K 

0.1652 

0.09345 

0.28916 


Table 3.1 : Universal KPZ values for cumulants of height distributions in d = 1 + 1 [28]. 



Flat 

Curved 

Flat-Growing 

ID Groove 

Stationary 

(X)c 

-0.75(5) 

-2.3(1) 

-2.4(2) 

-1.47(2) 

== 

(X 2 )c 

0.23(1) 

0.33(2) 

0.34(2) 

0.249(4) 

0.46(2) 

ISI 

0.423(7) 

0.33(1) 

0.33(2) 

0.396(7) 

0.244(8) 

K 

0.344(9) 

0.212(7) 

0 .21(2) 

0.31(2) 

0.176(4) 


Table 3.2: Universal KPZ values for cumulants of height distributions in d = 2 + 1 
dimensions for Flat case [57], Curved [58], flat-growing [60], 1D groove initial condition 
[58] and Stationary state [58]. 

As a final remark on KPZ height distributions, we point out that studying KPZ growth 
on enlarging flat substrates, Carrasco et al. have revealed that the Tracy-Widom distributions 
and the Airy processes (as well as their (2 + l)-dimensional analogs) do not depend on the 
interface macroscopic curvature, but actually on the inflation of the lattice metric on the active 
zone [60]. Furthemore, in the search for an upper critical dimension in KPZ class, Alves et 
al. [127] have confirmed that the KPZ ansatze is valid up to d = 6 + 1 dimensions, at least for 
the Kim-Kosterlitz model [128]. 
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3.2 Universal Squared Roughness Distributions 

In 1994, investigating the old problem of random-walk interfaces, Foltin et al. [129] showed 
that the squared roughness pdf [P(w; 2 )] 5 of such interface, at steady state , behaves as: 



where w 2 is the squared global roughness (eq. 2.9) of the interface, <b(tt) has a closed form 


and is also an universal scaling function. Indeed, its universality was confirmed by numerical 
simulations of EW and KPZ models, and have emphasized the power of that distribution for 
accessing the UC of a given growth process, once it presents a weak dependence on finite-size 
corrections. 

In the same year, Plischke et al. [130] extended theses studies for curvature-driven inter¬ 
faces and have reinforced the validity of the eq. 3.7 as well as the universality of <I>, which in this 
case has a different form from that for Gaussian interfaces. Following the same idea, Racz et 
al. [54] calculated numerically P{w 2 ) for several UCs in d = 2 + 1 dimensions and exposed that, 
unlike the lower-dimensional case, <&ew and &kpz are different, being the first a Gaussian and 
the last one marked by a slow decay in the rigth tail. <&mh and <&vlds also exhibited remark¬ 
able differences between each other. In the same paper one finds a “recipe” of how comparing 
4>. for a given UC, with the ones obtained in experiments. The recipe consistes in divide the 
surface into boxes of lateral size l <C £|| inside which w 2 should be calculated to yield a large 
ensemble. On the other hand, l must also be larger than characteristic sizes at surface such as 
grains, mounds, etc. 

The Squared Local Roughness Distributions (SRLDs) have been used in signal analysis 
[131] and in issues about the upper critical KPZ dimension [55] because they are considered one 
of the most suitable ways for accessing the UC of a growth. For instance, in a very interesting 
study [132], Aarao Reis analyzed the rescaled distribution 4>, at mean null and unitary variance, 
(call it T, defined as in the eq. 3.8) for one- and two-dimensional models belonging to the KPZ 
and VLDS classes. 



(3.8) 


Aarao Reis have confirmed that, ind = 2+ 1, ^ kpz presents a stretched exponential 
in the rigth tail as approximately exp(— x 0 ' 8 ) [132] - see figure 3.4(a). This form is contrasted 
with those from '&ew, which is Gaussian, and with the simple exponential decay of ’P vlds 
and ^mh- Even the last ones can be easily distinguished in the &vlds and <&mh scaling, as 
exemplified in the inset (B) of fig. 3.4(a). A subsequent work from Paiva & Aarao Reis [133] 
brought out '&kpz in the growth regime presenting a similar decay at the right tail as shown in 
the fig. 3.4(b). This feature has been suggested as an universal and distinct KPZ landmark [133] 
and has been confirmed experimentally in the growth of CdTe [63] and oligomer [64] thin films. 

5 Notice that at the steady state, the surface roughness fluctuates around its saturated value and P(w 2 ) is the pdf 

associated to these fluctuations. 
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Figure 3.4: (a) The main plot exhibits the contrasting form of universalV distributions 
for the VLDS and KPZ classes calculated from the CRSOS (L = 64), RSOS (L = 256) 
and RSOS2 (L = 128) models at the stationary state, for d = 2 + 1, according to the 
ref. [132], The insertion (A) reinforces the difference between 4/ vlds and ’f K pz near 
of the peak, whereas inserts (B) shows the constrast between §vlds and $mh- (b) 
Comparasion between 4 / KPZ calculated from the steady state (solid line) and from 
two-dimensional RSOS and BD models in the growth regime, for a “box” of lateral 
size r = 64. Results in (a) and in (b) can be found in the references [132] and [133], 
respectively. Theses results were kindly provided by Prof. Aarao Reis. 


3.3 Universal Maximal Relative Height Distributions 


Extreme-value statistics (EVS) play an important role in systems where rare events have drastic 
consequences such as floods, internet failures, stock market crashes as well as in important tech¬ 
nological applications [56, 125,134-137]. For instance, the onset of a breakdown of corroded 
surfaces is determined by their deepest (or weakest) point, whereas in batteries the highest point 
of the metal surface reaching the opposite metal surface is responsible by the beginning of a 
short-circuit [134]. A well-known probability function in the EVS context is the called Gum- 
bel’s first asymptote, being the distribution of the nth value among N independent (uncorrelated) 
random variables [125]. The Gumbel pdf, G(X ; m), of the variable X is defined as: 


G(X- m) 


m m b 

T f (m) 


exp[-m( 2 .Y + exp(—£x))], 


(3.9) 


where m is a parameter, b = \/ipi(rn)/(X 2 ) c , zx = b((X) — X + s), s = [ln(yn ) — 
r ipo(m)\/b. Vf(X) is the gamma function, and 'ipk(X) is the polygamma function or order k 
[56,125]. 

Early studies on EVS applied on growing interfaces focused on the steady state of lin¬ 
ear (EW) equations [134-136]. Through maximal relative height (rn*) analysis, defined as the 
difference between the largest height minus the average height of the surface, it was shown that 
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in the stationary regime m* scales likewise the global roughness (this result is valid only for the 
one-dimensional case) [134] and that its universal distribution, P(m *) = L -1 / 2 f(m*L~ 1 / 2 ), 
is dictated by the Airy distribution function, f(x), whether periodic boundary conditions are 
used [136]. 

Interestingly, even for strongly correlated systems such as two-dimensional EW inter¬ 
faces at steady state, it has been shown that P(m*) can be very well fitted by the Gumbel 6 , with 
a non-integer value, n = 2.6(2) in this case [137]. A numerical work has brought theses discus¬ 
sions for KPZ and VLDS two-dimensional surfaces [56], It was confirmed the universality for 
either the scaled maximal- or minimal-relative height distributions (MRHDs), depending on the 
relevant nonlinear term signal. Moreover, the remarkable contrast between the right tail decay 
from MRHD^pz (simple exponential) and MRHDyv.DS (Gaussian) and the very weak finite- 
size effects affecting these distributions have suggested an alternative way for accessing the UC 
of a given growth. Our experimental results have been the first confirmation of such universality 
of this distribution in the KPZ^ = 2 +i context. 


6 See also the Oliveira et al. results [59] about the relationship between the Gumbel distribution and KPZ 
HD’s<j = 2 +i. They are already discussed in the section 3.1. 




Material and Experimental Methods 


In this chapter a short review on the experimental techniques used along this work is given. 
Details on growth parameters and on the experimental methodology is discussed in details. 


4.1 Hot Wall Technique 

Hot Wall Epitaxy is a well-established technique based on thermal evaporation developed at the 
end of the 1970’s which has been used for growing high-quality films from II-VI, IV-VI and also 
III-V compounds [8,138]. The great difference between a simple thermal evaporation system 
and a HWE one consists the presence of a heated liner (hot wall). This liner serves as a guide 
for the vapor beam to flow from the source towards the substrate, ensuring low material loss and 
growth conditions as near as possible of the thermodynamic equilibrium. 

HWE systems share the same basic structure (see fig. 4.1), but they can be modified 
depending on the particular growth necessities (see [138] for different HWE forms). In general, 
there are tree resistances windings at quartz tube to heat the substrate, the source and the wall, 
independently. It guarantees that the temperature displayed at the controller is always being 
measured at the same referential 1 , beyond of making the system’s calibration to be very reliable. 
HWE works in high vacuums (« 10“' Torr), providing a relatively clean environment for the 
growth [8, 139, 140]. Substrate position can work as a lip for the wall or can be slightly inserted 
above it when a shutter is available. 

As a very simple technique, HWE is unsuitable for doping, for growing compounds with 
very different vapor pressures and it also does not allow in-situ measurements. Indeed, these 
disadvantages turn the HWE use minimized. However, in situations where in-situ measurements 
are not essential, where the compounds evaporate congruently and do not have high temperature 

1 It is not the case for other systems where temperature measurements must be routinely calibrated, as occurs in 
MBE systems. 
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of sublimation, HWE is absolutely one of the most suitable techniques to be employed due to 
its high reproducibility, wide growth rate range (0.01 - 10 A/s) [8] and its relative low cost 
(ownership and maintenance). This is, particularly, the case for the growth of CdTe [8,9, 15, 
139-141], 



Figure 4.1 : Hot Wall system of Epitaxy Laboratory, Physics Departament of the Uni- 
versidade Federal de Vigosa, Brazil. The legend in the figure shows two furnaces 
(substrate and source) that work independently and a shutter which allows start¬ 
ing/finishing the growth. The pressure during the growth can reach ~ 10~ 7 Torr. 
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4.2 Atomic Force Microscopy 

After the 1986 year, AFM invention [142] has enormously assisted the surface study in microm¬ 
eter scale beyond of helping several stages of the thin-film production [3]. As a blind man must 
use his fingers to feel topographic variations during a Braille lecture, the AFM consists in to 
obtain morphological data from interactions between a very thin tip (order of 10 _8 m radii) and 
the sample surface to yield a three-dimensional image. The equipment measures atomic forces 
via a spring deformation, more known as cantiveler, under which the tip is coupled. On the 
other cantiveler face, a laser beam is focused and reflected to a photo-diode in order to send 
electrical signals associated with the cantiveler’s deflection to the controller. A simplified view 
of most important AFM components is shown at fig. 4.2. 

Once the electrical signals have arrived at controller, they are translated and sent to the 
computer to yield a bidimensional profile sketched by the tip. An AFM image can have 1024 
x 1024 pixels, which means that there are 1024 bidimensional profiles, each one of them being 
composed by 1024 equally spaced points. In some AFM systems, the sample is supported by 
a piezoelectric ceramic which suffers well-behaved subnanometric deformations to move the 
sample in relation to the tip. 

There are two basic AFM measurement modes, namely, contact and tapping. In the 
former, the tip approaches to the surface until suffer a repulsion force, which causes the tip to 
bend up. In the tapping mode, however, the tip is set to oscillate near to its natural resonance 



Figure 4.2: Basic scheme of working principles of an Atomic Force Microscope. 
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frequency and gets close to the surface until the oscillation amplitude becomes reduced to a 
smaller referential value (this occurs at distance between regions of repulsion and attraction). 
For both modes, the piezo adjusts, through a feedback mechanism, the height axis (z) to hold or 
the force (contact) or the amplitude (tapping) constant. Then, a surface image is obtained from 
those z values recorded by the piezo. The choice of which method should be employed depends 
on the aspects of the studied material (shape, structure, nature) and the kind of results that 
want to be found (local friction coefficients, phases difference, conductance variations, etc). For 
instance, biological samples are most suitable for tapping mode once it avoids damage sample 
and frictional forces, while contact is most suitable for thin films and crystals in general. 

Nowadays, AFM is one of the most used techniques for studying surfaces at the submi¬ 
crometer level. In particular, several works have used AFM images to perform scaling analysis 
of interfaces, as exemplified by works involving the growth of Si02 by CVD [61], the disso¬ 
lution of polycrystalline pure iron [85], the amorphous Si by thermal evaporation [97] and Pt 
sputtered on glass [99]. 


4.3 Si surface cleaning 

Substrate cleaning is the key initial step in thin film and epitaxial growth. In order to obtain 
flat and contamination-free Si surfaces, several chemical and/or thermal treatments have been 
proposed since 1980’s [143-145]. Basically, the aim is to remove the silicon native oxide layer 
(« 0.7 nm of thickness) and the hydrocarbon contaminant layer (~ 0.2 nm). The former is 
usually removed by low-concentration aqueous F1F solution [1-10%] using F1F and water of 
high purity, which provide an “ideal” stable monolayer Fl-terminated Si surface, depending 
on the aqueous pH concentration [146-148]. Moreover, aqueous HF solution does not etch 
the bare Si surface itself, preserving a Si smooth morphology (typically, ~ 0.2 nm of rough¬ 
ness [145]). Whereas some works stress the necessity of removing reminiscent hydrocarbon 
impurities [143,144], X-ray photoelectron spectroscopy (XPS) measurements show that 1.5% 
HF-treated Si surfaces present a very low concentration of O, F and C [145]. Moreover, the 
H-terminated surface obtained by this process proved to be very stable against the oxidation in 
air. In particular, this cleaning procedure has been used for growing epitaxial CdTe QD’s on 
Si(l 11) [ 8 ], Another approaches as the degreasing by acetone, ethanol and deionized (DI) water 
followed by repeatedly boiling in HNO 3 , dipping in HF, rinsed with DI water and dried with N 2 
also are commonly used [139]. 
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4.4 CdTe thin films: Cleaning, Growth and Characterization 

In this work, p-type Si(001) substrates of dimensions 10.0 mm x 10.0 mm x 0.3 mm were 
dipped 2 minutes into an aqueous HF(2%) solution prepared with DI water. This time has been 
checked to be able of removing completely the Si native oxide layer [145], At the sequence, 
Si surfaces were exposed to N 2 air jet just to remove any reminiscent droplet on the surface, 
once after HF treatment the surface becomes largely H-terminated. The treated surfaces were 
immediately inserted into the HWE chamber, inside which a high vacuum was performed. There 
was no any heating treatment before or after the growth. 

The HWE system used in this work is composed by two independent furnaces (source 
and substrate), separated by a shutter and a wall of 7 cm, as shown in the fig. 4.1. The depo¬ 
sition occurs at pressures ~ 10 ' Ton - , obtained by a diffusion pump system, while the source 
temperature can be controlled from 400 to 520 °C, producing growth rates between (~ 0.01 and 

2.5 A/s). The temperature of the substrate (more known as the deposition temperature) can be 
varied from 150 up to 550 °C. 

Solid CdTe (99.999%) has been used as source material. The temperature of the source 
was fixed at 520 °C yielding a deposition rate F = 2.2 ± 0.3 A/s for deposition temperatures set 
at 150, 200, 250, and 300 °C. For each deposition temperature, the growth time (t) was varied 
from 15 to 240 min in a geometric progression sequence of ratio 2, providing thicknesses (th) 
from approximately 2 0.20 up to 3.5 pm. The th and the growth rate F were determined post¬ 
growth using a XP1 - AMBIOS contact profilometer and a ContourGT-K BRUKER optical 
profilometer. 

Surface characterization was performed in air by ex-situ AFM. We have used an Ntegra 
Prima SPM working in contact mode. Different kinds of Si tips were used in order to check the 
reliability of data. All of them were used in the statistical analysis. The frequency of the AFM 
scan was kept near of 1.5 lines/s during acquisition of the images, but we have also confirmed 
that the frequency does not affect the results as far the frequency is not set at very high values 
as > 4.0 lines/s. Surface topography of 3 to 10 different regions near of the center film of 
each sample was scanned producing images of 10 pm x 10 pm with 1024 x 1024 pixels. 
This size was chosen so that morphological properties in domains smaller and larger than the 
averaged grain size could be simultaneously investigated. However, different sizes of scanning 
were earned out, namely, 1 pm x 1 pm, 30 pm x 30 pm and 100 pm x 100 pm to guarantee 
the validity of data in different scales. For all images a flatten’s correction of second order was 
performed to correct the sample misalignment and the piezo scanner error. 


2 As we shall see in the chapters 5 and 6, the growth of CdTe in these conditions are non-conserved, this implies 
that the growth rate is not constant in time. 
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RESULTS: Uncovering the KPZ Universality in CdTe 
Thin Films 


The results contained in this chapter are related to films grown at a particular deposition tem¬ 
perature, namely, T = 250 °C. Here we develop a noveil procedure to distill the Universality 
Class (UC) of a given growth, which consists in to: i) perform a visual-AFM investigation of the 
interfaces fluctuations focusing on the dynamic of superficial structures whose length is char¬ 
acteristic; ii) relate it to a deep local scaling analysis; iii) add on the results coming from the 
global scaling and, finally, iv) to find and/or to confirm the UC of the growth by using universal 
distributions. As we will be seen below, through this scheme we were able to find the UC of 
CdTe thin films, in the meantime that one has been experimentally demonstrated the universality 
ofKPZ d=2+1 distributions. 

5.1 Semi-Quantitative Morphological Analysis 

At first, it is important to observe how the morphology evolves as function of the growth time in 
a semi-quantitative fashion. Figure 5.1 shows 10 X 10 /im AFM images for all available times. 
The images depicted are those which better represent typical behavior over all regions scanned. 

At initial times (figures 5.1 (a) and (b)), the surface is dominated by a large number of 
grains with well-defined “sharp” 1 shape, which present a very high aspect ratio (Q), defined 
as the ratio between the crystallite height and the characteristic length of its base (the average 
grain width - Q. Roughly, one observes ( < 0.5 /im for t = 15 and 30 min, whilst the vertical 
bar decreases from 80 nm to 60 nm. Consequently, il seems to decrease in this interval. At 
the subsequent time, one can see a surface composed by fewer, but larger grains (£ ~ 0.5 lira) 

1 Hereafter we consider that the grain has a sharp shape under a [nm] x [gm] scale, which is a typical scale for 
analyzing the grainy morphology [12,15,61]. 
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X (|i/7?) 

Figure 5.1 : AFM images [10 x 10 /xm; height scale in nanometers] for CdTe thin films 
grown at T = 250°C. (a) - (e), Interfaces grown by 15, 30, 60, 120 and 240 min, 
respectively, (f) Typical grain’s shape present at the interface for t = 60 min (solid 
lines) and t = 240 min (dashed lines). 
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surrounded by a crowd of smaller sharp grains. It is not difficult to see that ft has decreased once 
again. Remarkably, for the last two surfaces showed in the figures 5.1 (d) and (e), one notices 
the presence of grains with larger basis, in particular, ( < 0.8 prn, and having a smoother top, 
contrasting with those previous conical-like grains. Considering that the surface (d) still presents 
a background of small crystallites, the aspect ratio in the interval of 120 — 240 min keeps near 
to a constant value, once that background has grown up both in height and width. Another 
important and distinguishable feature differing surfaces (d)-(e) from (a)-(c) is the ubiquity of 
multi-peaked grains at the interface. The AFM images point out that coalescence processes 
have yielded these structures, at least for times larger than t = 60 min. In the fig. 5.1 (f) one 
compares typical grain’s profiles at the surface for t = 60 min and t = 240 min. It is clear that 
single-peaked-sharp-conical structures evolve for larger multi-peaked mounds with a flatter top 
and having a smaller aspect ratio. 

The origin os theses results can be understood as follows. Due to the large mismatch 
of lattice parameters between CdTe and Si (~ 20.0% [8]), the growth of CdTe on Si hardly 
proceeds epitaxially, except in particular conditions reported by Ferreira et al. onto Si(lll) 
substrates [8,9]. Otherwise, CdTe layers are usually polycrystalline presenting a strong [111] 
texture [13, 141, 149]. Indeed, we confirm this by performing a structural analysis of the 
CdTe/Si(001) films. Using a D8 - DISCOVER (BRUCKER) X-Ray Diffractometer in the 9—29 
coupled mode (radiation A cuK a = 0.154056 nm), the XRD spectra of films are built, as shown 
in the fig. 5.2. The presence of more than one crystallographic orientation is a sufficient con¬ 
dition to assign the polycrystalline nature of film. Nevertheless, as the deposition proceeds, the 
(111) peak becomes higher at the same time that the other ones are reduced. The probability of 
finding a crystallite [111] oriented after t (min) of growth, p( 111, t), is defined as: 

„(iim) . /yy ; 11 < 5 .i) 

l^hkl 1 \ t )hkl/A g _2e{Vhkl) 

where Ihkl.t, is the intensity of the peak (hkl) at the time t and A g _2o(9hki ) is the absorption 
factor for the 0 — 20 geometry dependent of the O^m angle [4,150]. 

In the inset of the fig. 5.2 one shows the probability of finding a (111) grain in the film. 
Thus, one can confirm a very strong [111] texture in CdTe layers grown on Si(001) substrates at 
T = 250°C where for early growth times (111) crystals have already composed almost 90% of 
the thick layer. 

In regard to initial regime of growth, previous studies have shown that the growth of 
CdTe on FI-terminated Si surfaces follows the Volmer-Weber (VW) growth mode, where three- 
dimensional islands nucleate directly onto the surface without an initial wetting layer [8, 139, 
140]. This result has been found in the CdTe growth on both (111) [8, 139] 2 and Si(001) 
faces [140], by using a slow and/or a variable molecular flux during the experiment. Due to 
the VW growth mode, conical-like grains are formed onto the template of submonolayer islands 
and grow by increasing their aspect ratio, since the direct contact with the substrate is unfa¬ 
vorable energetically. Parts of the Si remains exposed until the islands start coalescing to form 
a continuous layer. As the environment is polycrystalline, grain boundaries (GB) of colided 
neighborhing grains are unfavorable regions for deposition and diffusion of particles, since in 

2 See fig. C.4 in the appendix section C.2.1. 
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Figure 5.2: X-Ray spectra of CdTe polycrystalline layers using \c U K-a ■ Spectrum for 
layers grown by (black) 15 min, (red) 120 min, and (blue) 240 min at T = 250 °C. One 
notices the presence of grains with several crystallographic orientations in reference 
to the substrate normal, namely: (111), (220), (311), (400) and (331). However, as 
the time evolves, the (111) peak becomes higher. The inset shows the probability of 
finding a grain whose [111] direction coincides with the normal vector of the substrate 
surface. 


these places one finds a large number of defects [151-154]. Defects are yielded, mainly, due 
the stress evolution. Many studies [151,152] have shown that the stress is compressive prior 
to coalescence, becomes tensile during the coalescence (as a consequence of attractive forces 
between neighboring grains) and reaches a constant value as deposition proceeds. Even in the 
presence of the strong texture, the shared GBs of two (111) crystals also yield defects because 
the grains have usually different rotational orientations. This kind of defect causes shear stress 
(torque) which tends to reorient the crystallites [152]. 

For long growth times (multi-layer regime), the substrate-film interaction is expected to 
be lost [155] and the growth dynamic becomes dictated by the growth of CdTe on CdTe, i.e. 
the homoepitaxial growth [155, 156]. As the substrate temperature is high during the growth, 
the rate diffusion is expected to be high [22] 3 . Thus, “defect sites” at (an around) the GBs of 
colided neighboring grains are eventually covered by diffusion and deposition of particles (a 
kind of relaxation process occuring at, and around, these regions), which allows coalescence 
processes becoming more operative. As consequence of this relaxation, the aspect ratio (fl) of 
grains decrease in time and the top of them becomes smoother (due to the “filling” process of 
particles at local minima), as corroborated by the AFM images (figures 5.1(a)-(e)). 


3 See also section C.2.1 in the Appendix. 

















5. Uncovering the KPZ universality in CdTe thin films 


38 


5.2 Local Fluctuations: Non-Universal and Universal Scaling Expo¬ 
nents 

Now we turn our attention to local fluctuations at the interface. Figure 5.3(a) shows the local 
roughness (defined in the section 2.2) calculated for surfaces grown by different times. 



Figure 5-3i ( 3 ) Local toughness scaling for CdTe thin films grown at T — 250 °C by 15 
min (left triangles), 30 min (up triangles), 60 min (diamonds), 120 min (squares) and 
240 min (circles). Inset shows the average local inclination and (b) shows in details 
the non-universal (geometrical) scaling exponent as function of the growth time. Solid 
lines in the main plot of (a) represent the fits used to extract on in (b). 


At short length scales (l < 10” 1 //to) and early times, one can observe wi oc (l. t) de¬ 
creasing in time because, in this situation, wi oc measures the intra-grain roughness of grains 
evolving from a shaip to a rounded shape - i.e, from larger to smaller height fluctuations. For 
comparison, this is the opposite occurring in systems that exhibits anomalous scaling, as can be 
seen in the appendix section B. In general, the presence or absence of anomalous roughening 
can be evaluate by calculating the evolution of squared slopes at the surface, ((V/ 1 ) 2 ), which is 
expected to scale as [69] : 


((V/ 1 ) 2 ) ~ t 2K . (5.2) 

Positive values of k indicate anomalous scaling, while negative or null values confirm 
that the Family-Vicsek scaling (eq. 2.11) describes the dynamic of local height fluctuations 
[69,70], 

The inset of the fig. 5.3(a) shows the behavior of ((V/ 1 ) 2 ) in time, which is the same 
of wio C for fixed box sizes /*, with l* < 10” 1 /im. In particular, one can notice that an initial 
decreasing regime is followed by a “saturated” one. The power-law (eq. 5.2) in this transient 
regime provides k = —0.15(5). This value is very close to that predicted by Lopez [69], namely 
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k = —1/6, for two-dimensional KPZ systems. Moreover, it is well-known that wi oc (l*,t ) ~ 
[157,158] and these results are in agreement with our previous analysis based on AFM images. 
Moreover, the crossover in wi oc occurs in l « ( (compare with tig. 5.1), in accordance with 
numerical results [76,77]. At the other hand, long-wavelength fluctuations (Z 3> 10 -6 pm) 
increase in time, as expected [22], 

The scaling exponent a\, defined as w(l, t) ~ V" , has its value changing in time from 
0.6 to 0.9, respectively, as shown in the fig. 5.3(b). These values, however, should not be taken 
as representative of universal fluctuations at the interface, as done in previous studies [12, 14,15], 
Rather than, as explained by Oliveira and Reis [76,77] (see section 2.3.2), a\ values are related 
to the grain morphology: a± ~ 0.6 indicates that the grainy morphology has a sharp-conical 
form and, as larger is the aq value, smoother is the top of typical grains’ shape. Notice that 
these findings are in agreement with the AFM images (fig. 5.1). The true universal roughness 
exponent a, however, should be found in a region where ( <C l <C £ [76,77]. 



Figure 5.4: Normalized slope-slope correlation function for CdTe thin films grown at 
T = 250°C by 15 min (left triangles), 30 min (up triangles), 60 min (diamonds), 120 
min (squares) and 240 min (circles). Inset shows the first minimum (circles) extracted 
from T(Z)/T(0) curves as function of growth time. 


In order to unearth the dynamic exponent, we have used the slope-slope correlation func¬ 
tion, defined in the eq. 2.15 and discussed in the section 2.2. Figure 5.4 shows normalized T(Z) 
functions calculated for all times studied. From the first zero and/or from the first minimum 
(r m ) of normalized T(Z) curves one can measure the coarsening exponent n coar , defined as 
r m ~ t ncoar [15, 63, 85]. For short growth times, before the appearence of multi-peak struc¬ 
tures at surface, one has that n coar = 1/z. As suggested by AFM images, this regime in our 
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experimental situation occurs for t < 60 min and provide n coar = 1/z = 0.62(2) - see the 
dashed line in the inset of the fig. 5.4. This value is in excellent agreement with that expected 
for two-dimensional KPZ systems [50,51,57-59,64] (see eq. 3.5). In turn, taking longest avail¬ 
able growth times, one obtains n coar ~ 0.34, which should not be interpreted as representative 
of temporal universal fluctuations because, here, peaks of reasonable size appearing on the top 
of a mound induce an underestimative for the average grain size measured by the slope-slope 
function (eq. 2.15). Thus, a reliable measurement for r m is not achieved in this situation and, 
hence, one can not be validated the relation n coar = 1/z. 


5.3 Global Scaling 

Following the Family-Vicsek ansatze (eq. 2.11), the global roughness scale as t :i . In the fig. 
5.5 one can see the w plot behaving as a power-law with exponent ft = 0.24(4). This is a very 
strong evidence of KPZd= 2 +i growth, once ftxpz ~ 0.24 (see eq. 3.5). 

At this point we have found 1/z = 0.62(2), ft = 0.24(4) and k = —0.15(5). These 
values are strong indicatives of two-dimensional KPZ growth in CdTe films. If this clue is really 
true, we should obtain universal distributions of heights, local roughness, and maximal heights 
matching with that ones numerically predicted for the KPZ class (see table 3.2 and sections 3.2 
and 3.3). 



Figure 5.5: Global roughness as function of time for CdTe thin films grown at T = 
250 °C. The roughness grows as a power-law with ft = 0.24(4). 
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As we arc dealing with a few number of available growth times, we are not able to find 
the asymptotic velocity (Voo) of the CdTe surface growth [113]. Beyond that, by following 
the Krug-Meakin toolbox [113,118], data of CdTe interface fluctuations evolving in time (until 
reaching the stationary regime) for different initial miscut angles of Si(001) substrates should be 
accessible. Thus, we are not in position to validate, at this time, the two-dimensional KPZ an- 
sdtze (eq. 3.6). However, we can perform a comparison between a rescaled height distributions 
and the distribution of x- Take the spatial average on both sides of the eq. 3.4 - negleting the 
correction terms - one has: 


(h) = Voot + sign(\)(Tt) 0 (x) ■ (5.3) 

One can vanish the term t subtracting eq. 5.3 from eq. 3.4. The result is: 


h-(h) = sign(X)(Tt)P[x - (*>]■ (5.4) 

Now, it easy to show that (h 2 ) c = (Tt) 2 / 3 [(x 2 ) — (x) 2 ]- Hence, dividing eq. 5.4 by the 


global roughness w = Oh = ( h 2 ) l J 2 one finds: 


(h-(h)) 


sign(X) 


(x- (X)) 

a x 


(5.5) 


where a x = (x 2 )c /2 - 

Equation 5.5 means that to compare rescaled height distributions with mean null and 
unitary variance is equivalent to comparing the x fluctuations rescaled in the same way. This 
allows us to compare these distributions even without known the non-universal KPZ parameters. 
Figure 5.6 depicts the rescaled height distributions for several regions of the thickest sample. 
Remarkably, one notices a very good agreement between the experimental data and the Gumbel 
pdf with m = 6 (see eq. 3.9) for four decades around the peak. This astonish result is reinforced 
by the data collapse close to the peak, as highlighted in the the insertion A of the fig. 5.6. The 
slight deviation observed in the left tail could be due to the fact that AFM tip does not scan 
accurately deep valleys as well as it does for higher heights. 

Comparing the universal KPZ f / = 2 +i dimensionless cumulant ratios for the flat case (table 
3.2), we have found, for the thickest film, St= 24 o = 0.34(1) and K t= 240 = 0.3(1), while 
numerical values are S = 0.423(7) and K = 0.344(9) [57-59,64]. However, as shown in the 
inset B of the fig. 5.6, the experimental S and K values approach to the KPZ,/ =2+ i ones only 
for t = 240 min. In principle, these results could be just a fluctuation, making the agreement 
of HDs in the fig. 5.6 a mere coincidence. Indeed, this is not the case here because before to 
seek the universal HD distribution, critical exponents have already point out to KPZ growth. 
Moreover, it is well know that experimental S and K values take a long time to converge to 
the asymptotic values calculated numerically [38,39,46]. This effect should be stronger in 2 + 
1 dimensional systems. Hence, the underestimate S value can be only a finite-time effect and 
results emerging from the fig. 5.6 provide the first so-waited experimental confirmation of the 
KPZrf = 2 +i universality beyond the exponents. 

In order to prove the universality of distributions described in the sections 3.2 and 3.3, 
and also to gain further evidence of KPZ growth, we have calculated the SFR and MRH distri¬ 
butions, constrained to the box size l within the interval [10/im/1024] < I C (. Figure 5.7 (a) 
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[h - <h>]/a h 

Figure 5.6: Rescaled height distributions at mean null and unitary variance from sev¬ 
eral regions (symbols) of the thickest sample grown at T = 250 °C compared with the 
numerical KPZ curve (solid magenta line). Inset (A) highlights the very good collapse 
between experimental and numerical data close to the peak in a linear x log plot. 
Insertion (B) exhibits the experimental S and K evolving in time. Solid and dashed 
magent lines refers to the numerically expected KPZ d=2+ i S and K values, respec¬ 
tively. 


shows the SLRDs for all available growth times, whereas in fig. 5.7 (b), the respective S and 
K values as function of l are depicted. Astonishingly, a very good collapse is seen, in almost 
four decades around the peak, among the experimental and numerical KPZ SRLDs. Moreover, 
the stretched exponential decay at the right tail is present, giving one more (strong) evidence of 
KPZ growth. 

Similar results are found for the MRHDs, figures 5.8 (a) and (b), where the nice collapse 
between the experimental and numerical data is evident. The curves showed in the main plot are 
those whose S and I\ values better represent the mean for the l described in the graph. 

It is very important mention why the S and K values change with /. This happens 
because one does not know, a priori, what are the Vs satisfying the condition ( < ! « ( 
(remember there is a characteristic length at the surface). As in the CdTe/Si(100) system ( « £, 
the appropriate interval for l is very short and should be identified as a convergence of S and 
K for some specific value. In the case shown in the figures 5.7(b) and 5.8(b), these convergent 
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Figure 5.7: (a) Rescaled squared local roughness distributions for samples grown at 
T = 250 °C, by 15 min (orange left triangles), 30 min (blue up triangles), 60 min (green 
diamonds) and 240 min (black circles).The curves shown are those whose the box size 
is indicated below the legend. They are different for each t because the appropriate 
interval C < l < £ depends on t. The inset shows the good collapse around the peak 
for the thinnest and thickest film grown, (b) Skewness and kurtosis values as function 
of the box size l. 
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Figure 5.8: (a) Rescaled maximal relative height distributions for samples grown at 
T = 250°C, by 15 min (orange left triangles), 30 min (blue up triangles), 60 min 
(green diamonds) and 240 min (black circles). The curves shown are those whose 
the box size is indicated below the legend. They are different for each t because the 
appropriate interval C < l < £ depends on t. The inset shows the good collapse 
around the peak for the thinnest and thickest film grown, (b) Skewness and kurtosis 
values as function of the box size l. 
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values agree with those numerically calculated for the SLRD and MRDS of KPZ,/ = 2+1 models 
[56,132,133]. 

5.4 The origin of the KPZ mechanism and Conclusions 

The results presented in the previous sections point out that the long-wavelength CdTe surface 
fluctuations, at least for film grown at T = 250 °C and F ~ 2.2 A/s, evolve according to the 
KPZ equation (eq. 3.1). The origin of the KPZ scaling in this experimental system can be 
understood unveiling the complex dynamic of grains, which is dictated by an interplay between 
inter-grain surface energetics, constraints yielded by the coalescence of neighborhing grains and 
filling processes occuring at (and around) the GB’s of colided neighborhing grains. This leads to 
a complex packing of crystalline grains. A simple illustration is provided by the grain deposition 
model in the Ref. [76]. In the model, a new cubic grain is firmly attached to the boundary 
of the grains below it, but does not fill all available space in their neighborhood because its 
shape is constrained to be cubic. This aggregation mechanism has the same effect of the lateral 
aggregation as in the ballistic model [110]. It generates excess velocity, which is the landmark 
of KPZ scaling. 

In conclusion, we have found the first robust 4 experimental KPZd= 2 +i system. On 
the basis of a semi-quantitative morphological and long-wavelength fluctuation studies, we have 
extracted critical exponents agreeing with that expected for the KPZ class. Moreover, the univer¬ 
sality of KPZ distributions (HDs, SLRDs and MRHDs) has been experimentally demonstrated. 
We have not checked the KPZ c i= 2 +i ansatze (eq. 3.4) directly due to inherent two-dimensional 
experimental obstacles. Our results were published in the Physical Review B [63], 

Shortly after the publication of our work [63], T. Halpin-Healy and G. Palasantzas [64] 
have used this noveil scheme proposed by us to confirm the KPZ universality in oligomer thin 
films. This clearly corroborates that we have given a new perspective for finding, and confirm¬ 
ing, two-dimensional growing surfaces belonging to the KPZ universality class [64]. 


4 Going beyond of the standard comparasion with exponents. 





RESULTS: The Effect of Temperature on CdTe Growth 
Dynamic 


In this chapter one studies the effect of the deposition temperature, in the range of [150, 300] °C, 
on the mound evolution and on the long-wavelength fluctuations of CdTe surfaces grown on 
Si(001) substrates. 


6.1 Intragrain morphology and local fluctuations 

Figure 6.1 shows typical AFM images for surfaces grown at different T for the two largest 
available growth time, t = 120 min and 240 min. For T = 150 °C (fig. 6.1(a) and (b)), one 
can see grains with a well-defined sharp shape dominating the surface, at the same time that Q 
seems to increase with t, as suggested by the vertical bar. For films grown at T = 200 °C, the 
scenario is quite si mi lar to that for T = 150 °C, except for the largest growth time (fig. 6.1(d)), 
where has decreased in the interval of t £ [120, 240]min. Unlike the morphology of surfaces 
grown at T = 250 °C (fig. 5.1), multi-peaked structures arc rare in t = 240 min and grains with 
a well-defined shape still are the majority. See fig. 6.2(a), in which typical mound profiles arc 
shown in order to make clear those notes. On the other hand, at higher temperatures, namely 
T = 300 °C, the surfaces present much more complex structures, as can be noticed from the the 
figures 6.1(e) and (f). In fact, these films arc composed by a mix of a few conical grains and large 
mounds (formed by coalesced grains) with £ ~ 2 fan. As the time evolves, Q increases faster 
than the mound height, at least within the interval t £ [120, 240]min, to form larger mounds. 

The orientation of grains whose crystallographic plans arc parallel to the substrate sur¬ 
face is determined by 0 — 20 XRD measurements. As we already found for films grown at 
T = 250°C, a strong [111] texture has been observed also for other temperatures. Moreover, 
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Figure 6.1 : AFM images (10 x 10 nm, height scale in nm) of CdTe thin films grown at: 

T = 150°C (a) and (b), T = 200°C (c) and (d), and T = 300°C (e) and (f) by t = 120 
min and 240 min, respectively. 

the same peaks have been seen in the XRD espectra 1 . Figure 6.3 shows the probability of find¬ 
ing (111) grains (eq. 5.1) in the CdTe layer as function of the growth time and for different 
temperatures. The results suggest that this texture is: 1) slightly influenced by the deposition 
temperature, being this sligth deviation probably due to the fact that as higher is T, larger is the 
width of grains [159], including those grown in (11 l)-different orientions which reflects more 

'The XRD spectrum is not shown because it is similar to that represented in the fig. 5.2. 
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T = 300 °C for different growth times. 



Figure 6.3: Probability of finding (111) grains with respect to the substrate normal in 
a CdTe layer grown at T = 150°C (circles), 200°C (squares) and 300 °C (diamonds). 


the X-ray incident and, hence, diminhes p(lll,t,T). 2) independent of the substrate, once a 
similar [111] texture has been found by Ribeiro et al. [13] and Ferreira et. al [8, 141] in the 
growth of CdTe on glass substrates and also by Sporken et al. using Si(001) substrates [149]. 
Anyway, the important for the growth dynamic is that most coalescence process happens be¬ 
tween (111) grains and, even this way, they give rise a large number of defects at the GBs of 
colided neighboring grains due to, basically, the difference between the rotational orientation of 
neighboring crystallites. 

In regard to surface fluctuations at local scales, fig. 6.4 presents the local roughness 
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Figure 6.4: Local roughness for films grown at (a) T = 150°C, (b) 200°C and (c) 
300°C by 15 min (left triangles), 30 min (up triangles), 60 min (diamonds), 120 min 
(squares) and 240 min (circles). Insertion of (b) shows the squared local inclinations 
as function of the growth time for T = 150°C (circles), T = 200°C (squares), T = 
300°C (diamonds). The same legend is used in the inset of the figure (c), which 
depicts the geometrical a x exponent as function of t, extracted from the solid lines 
shown in the main plots. 


calculated as function of t and T. The first note to take look at is that, in the fig. 6.4(a), the 
curves dislocate to up as the time evolves (in both short-, l < 10“ 1 /im, and large-wavelengths 
l 3> If) " 1 firri) and, apparently, no one relaxation is observed. Indeed, this wi oc increasing 
at short-length scales is the landmark of the anomalous roughening, discussed in the appendix 
section B, where, at first sight, a larger number of critical exponents must to be found for deter¬ 
mining the UC. The k exponent, defined in the equations 5.2 and B.9, has been calculated for the 
films grown at T = 150 °C - see insertion of the fig. 6.4(b). The value obtained is k = 0.5(1). 
Although this value is positive, indicating anomalous scaling, its particular value is consistent 
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with that one obtained from an uncorrelated interface, where the squared-local-slope fluctua¬ 
tions increase in time with unit power. In another words, the random growth is intrinsically 
anomalous in this sense. 

For surfaces grown at T = 200 °C, the situation is more complex than that for the lowest 
T, as shown in the fig. 6.4(b). There, one can see that the curve for t = 240 min has dislocated 
down at short-length scales, instead of dislocating up, as in the initial trend. Very interesting, 
this has been the same behavior presented by fi, as well as in the visual inspection of AFM 
images (fig. 6.1). The squared local slopes also show a behavior of two regimes (see inset of 
the fig. 6.4(b)), where an initial anomalous regime, characterized by k = 0.19(5), is followed 
by a second one in which the Family-Vicsek scaling is asymptotically recovered [69, 82] with 
re » -0.7. 

For the deposition temperature T = 300 °C, the same qualitative result of two regimes 
is observed - see fig. 6.4(c). The transient anomalous scaling gives place to a “normal” one at 
f < 60 min, as corroborated by the local roughness at short-length scales and by ((V/;) 2 ) as 
function of t. The values for k characterizing these two regimes arc, respectively, 0.2(2) and 
—0.56(3). Notice that the first value can not guarantee the presence of anomaly, since k = 0 is 
a possible value within the error bar. 

Now, regarding to the geometrical scaling exponent a\, defined as wi oc ~ 1°' (section 
2.3.2), its value as function of the growth time is depicted in the inset of the fig. 6.4(c) for 
different temperatures. a\ is close to 0.6 for 150 °C, which indicates the presence of very sharp 
grains at surface at all times [76,77], according with the visual inspection of AFM images. In 
turn, for T = 200 °C and 300 °C, a\ changes from 0.6 to 0.8 and from 0.80 to 0.95, respectively. 
This reveals that sharp grains give place to smoother structures at long times. These results 
are totally consistent with the figures 6.2(a) and (b). A similar behavior for a\ occurs for a 
fixed t, as T increases. Based on these results, one can also conclude that the behavior of 
wi oc at short-length scales (/ <C 10 _1 jirn) is governed by two contributions: the first is the 
grain shape, related to a\, which can present large (sharp grains) or small (flat top grains) 
fluctuations. The second one corns from the aspect ratio, since for a fized box size l*, one 
has wi oc (l*, t) ~ ((V/i) 2 ) [157,158]. 

The origin of the results above can be understood as follows: initially, CdTe films evolve 
according to the Volmer-Weber growth mode [8, 139, 140], and then - width becomes larger as 
higher is the deposition temperature [159] (see fig. 6.1). As the grains with well-defined shape 
enlarge laterally, they colide forming a continuous film and giving origin to grains boundaries 
(GBs), where non-crystalline (defects) regions emerge [153,154]. These GBs framed by colided 
neighboring grains are formed mainly at the interface between (111) grains having different 
rotational orientations. The size of these defect regions is larger as higher is T since large are 
the grain perimeters (remember that initally £ increases with T) [159]. At these GBs, hence, 
an additional energy barrier ( Eqb ) rises, which tends to repel the diffusion and the deposition 
of particles at, and near to, these sites [151, 152]. A similar barrier has also been suggested 
recently in the growth of CdTe/CdS films [154]. For low T, where surface diffusion is slow, a 
small number of molecules can overcome the Eqb barrier and most of them aggregates inside 
the grain that they have arrived. This compels the grain height to increase faster than its width 
leading and ((V/i) 2 ) to increase as the time evolves (see figs. 6.1, 6.4 and inset of fig. 6.4(b)). 



6. Temperature Effect on CdTe Growth Dynamic 


51 


However, at long times, a relaxation happens because particles are eventually deposited at (and 
around) the GBs of colided neighboring grains. This diminishes the number of defect sites in 
the active zone and, consequently, enhances diffusion towards these regions. From this moment, 
coalescence processes become more operative and grains with well-defined shape give place to 
large mounds, so that <> and ((V/v) 2 } decrease in time. In particular, for T = 150 °C this regime 
is not observed due to the experimental condition of finite-time growth. For higher T, however, 
the second regime is clearly noticed and appears so earlier as higher is T, as can be checked in 
insets of the figs. 6.4(b) and (c) and in the inset of tig. 5.3(a). 

6.1.1 A Kinetic Monte Carlo Model 

The reliability of the above reasoning is illustrated in an one-dimensional atomistic growth 
model. Since our interest is the coalescence process, the growth starts on a periodic array of 
pyramidal grains with the same width ( and height H, for simplicity. The unit time is set to cor¬ 
respond to the deposition of one monolayer of particles. During a deposition event, a site i is ran¬ 
domly selected and a particle is allowed diffusing at surface (starting from the site i) until find¬ 
ing a site j where it aggregates permanently whether the constraint \hj — hj±i\ < 1 is satisfied. 
Thus, inside the characteristic “grain”, aggregation follows the conservative RSOS (restricted 
solid-on-solid) rule [160]. At the GBs of colided neighboring grains, however, an energy barrier 
Egb is present, so that a particle diffuses toward them with probability Pjj = e E cn/ k n T ^ 
Once a particle aggregates at a given GB i, the barrier Eq b at i becomes null with probability 
Pp> = fr E it/kn T \ Notice that this simple mechanism captures the underlying feature of the 
relaxation process at the GBs of colided neighboring grains. 



Figure 6.5: Results from an one-dimensional Kinetic Monte Carlo model. In (a) one 
depicts the height profiles for T = 150°C (top) and 200°C (bottom) for t = 10,100 
and 1000, and shifted by m = 10,80 and 960, respectively. Dashed lines represent 
the initial (t = 0) GB’s of colided neighboring grains, (b) Average-squared-local slope 
((V/i) 2 ) versus time for surfaces grown at T = 150°C (circles), T = 200°C (squares) 
and T = 300°C (diamonds). 
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Figure 6.5(a) shows typical surface evolutions for T = 150 and 200 °C. One has con¬ 
sidered Eqb = 0.10 eV, Eji = 0.30 eV, Q = 64 and H = 8, 16 and 24 for T = 150, 200 
and 300 °C, respectively. For T = 150 °C one observes grains with almost fixed width and 
increasing height. A similar behavior is observed at short times for T = 200 °C but, for large 
t, large mounds (formed by coalesced grains) appears. The same occurs for larger T. This 
qualitative agreement with the experiment is corroborated by the evolution of the squared local 
slopes displayed in the fig. 6.5(b). Comparing these results with the experimental ones (inset 
of fig. 6.4(c)), one can confirm that the inteiplay between the relaxation process at the GBs and 
initial conditions (Volmer-Weber growth mode and width larger for higher T), in fact, explains 
the CdTe/Si(001) mound evolution. Nevertheless, it is important mention that this simple model 
does not reproduce (and do not have this intention) the complex dynamic taking place in the 
growth of CdTe films, but it gives valuable insights about underlying features which rule the 
coalescence process. 


6.2 Universal exponents 

Now, once we have got an explanation for what is going on at local scales when the deposition 
temperature is changed, we are prepared to turn our attention to the analysis of long-wavelength 
fluctuations and to the Universality Class (UC) of the growth. 

Figure 6.6(a) shows the global roughness as function of the growth time for different 
temperatures. At first, one notices that w(T) increases with T. In fact, as the CdTe layer is 
polycrystalline, the complex competition between grains gives rise to inter-grain fluctuations 
that are larger as higher is the temperature because: 1) the grains, themselves, are larger and 
2) the valleys separating these structures are deeper as higher is T, as observed in the fig. 6.1. 
Similar results have been found in the growth of CdTe grown on glass substrates covered by 
fluorine doped with tin oxide [12,15]. 

Regarding to the growth exponent, for T = 150 °C, one obtains 3 = 0.51(4) (see 
inset of fig. 6.6(a)), which in association with k = 0.5(1) give us a strong evidence that one 
has a Poissonian growth at this low temperature. For T = 200 °C, however, one can notice 
that for initial times (t < 60 min), global roughness increases with initial ~ 0.5, but this 
regime is followed by a distinct one where, clearly, one finds a tendency of (5 becoming smaller 
- compare w(T) for T = 150°C and 200°C. Despite this fact, one finds (3 = 0.41(5) for 
T = 200 °C. It is worth mention that this value does not match with anyone known UC [22,25]. 
In particular this value has been found in homoepitaxial growth of Cu(001) at 250K [161], of 
Ag(lll) and Ag(001) at 300 K [162], and also in the growth of CdTe on glass substrates by 
HWEatT = 250 °C [12,15]. 

Finally, for films grown at T = 300 °C, experimental results point out 3 = 0.21(5), 
a value close to the KPZ one (eq. 3.5) [51, 57, 58] and also to the MH and VLDS classes - 
see section 2.3.1 and appendix section A.3, respectively. Indeed, diffusion should dominate 
the growth at very higher T, overcoming other mechanisms. Flowever, with only this evidence 
we can not determine which growth equation describes the surface fluctuations of the CdTe 
interface. 
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Figure 6.6: (a) Global roughness as function of the growth time for different tempera¬ 
tures, namely, T = 150 °C (circles), T = 200 °C (squares) and T = 300 °C (diamonds). 
Solid lines are guide to eyes, from which the (3 exponent has been found. Inset of (a) 
shows how the growth exponent depends on T. The solid line marks the value 0.239, 
which is consistent with KPZ d=2+ 1 [51,57,58]. (b) Slope-Slope correlation function 
for CdTe thin films grown at T = 200 °C. The oscillatory behavior is observed for all 
temperatures and the curve shown is a typical one, valid for all temperatures analyzed. 
Insertion in (b) depicts the first zero (r m ) of the Slope-Slope correlation function as a 
function of time for several T's. The legend is the same exhibited in the inset of (a). 
Dashed and solid lines are guide to the eyes, from where the n coar and 1/z exponents 
have been determined, explained in the text. 
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Figure 6.6(b) shows the slope-slope correlation function (eq. 2.15) for CdTe thin films 
grown at T = 200 °C. The same oscillatory behavior is found for samples grown at T = 150 °C 
and 300 °C (not shown here). As discussed in the section 2.2, from the first minimum or from 
first zero of T(l, t), one can extract the average mound size (r m ), which scales as r m ~ t Ucoar 
[15, 17,63,85]. In an appropriate range of time, one has n CO ar = 1 /z. The insertion in the fig. 
6.6(b) shows the first zero (r m ) as function of the growth time, for different temperatures. For 
T = 150 °C, multi-peaked structures do not appear at the surface and the equality n coa r = l/z 
must be valid for all t. One finds l/z ~ 0.07, which gives one more prove that the growth 
is uncorrelated at this T [22]. For T = 200 °C, two regimes can be reasonably seen, they are 
indentified as dashed and solid lines in the inset of fig. 6.6(b). From the first regime, one has 
l/z ~ —0.02 and from the second one obtains \/z = 0.6(1). The first value is consistent 
with random growth [22], whilst the second one agree with KPZd= 2+1 within the error bar (eq. 
3.5) [5 1,57,58]. Finally, for T = 300 °C, two short and long-time regimes are also observed and 
provide n coar = 0.32(5) and n coar = 0.7(1), respectively. The first regime is consistent with 
the l/z value for the VLDS class [22] and the second is near of the expected l/z value for the 
KPZ class. Both values agree with numerical ones whithin the error bar. Nevertheless, at this 
temperature we must be careful in assume n coar = l/z, once multi-peaked structures appear 
since initial growth times. For times t < 60 min, local peaks on the top of the mounds should 
contribute to underestimate the average size of the structure. However, at large growth times, 
local peaks become shorter compared with the basis of very large mounds at the surface (see fig. 
6.2(b)) and, hence, they should not strongly influence in the measurement of r m . Thus, the most 
careful action should be set n coar = l/z only for long growth times, at least for films grown at 
high deposition temperatures, as T = 300 °C. 


6.3 Partial considerations 

At this point, we have performed a scaling analysis of CdTe surface fluctuations based also 
in local features of the growth. From local roughness curves we detect transient anomalous 
scaling and crossover effects during the dynamic, which have been related to the emergence of an 
energy bander at the GB’s of colided neighboring grains and to the relaxation process. A Kinetic 
Monte Carlo model has supported our reasonings, reproducing qualitatively the experimental 
results. From the local roughness curves it was not possible to unearth the a exponent because 
the second (universal) regime has not been observed. However, the values found for a\ are 
consistent with the predictions in the refs. [76,77] as corroborated by the AFM images - see 
table 6.1. In the sequence, one has found /3 and l/z as function of the deposition temperature. 
For T = 150 °C, all results indicate a Poisson growth, while for T = 200 °C a crossover from 
uncorrelated-to-KPZ growth seems to occur. The value 3 = 0.41(5), which does not match with 
any known UC reinforces this idea. For the highest temperature studied, one has (3 = 0.21(5) 
and l/z = 0.7(1). The (3 value is close to the MH, VLDS and KPZ classes [22] and prevents 
a clear distinction of which growth equation describes the CdTe surface fluctuations at this 
temperature. However, the value found for l/z at long growth times points to a KPZ growth, 
even at this very high deposition temperature. Moreover, we stress that based only on results 
coming from exponents, one can not make clear what UC the fluctuations of CdTe surfaces (at 
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T = 200 °C and 300 °C) belong. 

The values of all exponents found until here can be seen in the tables 6.1, 6.2 and 6.3. 


t (min) 

15 

30 

60 

120 

240 

ai(T= 150°C) 

0.54(3) 

0.48(2) 

0.53(3) 

0.62(3) 

0.59(3) 

ai(T = 200 °C) 

0.61(2) 

0.58(3) 

0.63(2) 

0.74(2) 

0.82(2) 

ai(T = 250 °C) 

0.64(2) 

0.76(1) 

0.84(1) 

0.89(1) 

0.90(1) 

cri(T = 300 °C) 

0.78(1) 

0.87(2) 

0.89(2) 

0.92(1) 

0.93(4) 


Table 6.1 : Values for the non-universal exponent on (t, T ) coming from CdTe surfaces 
grown on Si(001) substrates by HWE with F « 2.2 A/s and deposition temperatures 

at 150°C, 200°C, 250°C and 300°C. 


T(°C) 

150 

200 

250 

300 

K 

0.5(1) 

0.19(5)/ « -0.7 

-0.15(5)/ « -0.05 

0.2(2)/-0.56(3) 

coar 

» 0.07 

-0.02(5)/0.6(l) 

0.62(2)/0.35(5) 

0.32(5)/0.7(l) 


Table 6.2: Values for the non-universal exponents «(t, T) and n coar (t, T) coming from 
CdTe surfaces grown on Si(001) substrates by HWE with F « 2.2 A/s and deposition 
temperatures at 150°C, 200°C, 250°C and 300°C. Values separated by the symbol 
(/) are valid in a appropriate range of time. These ranges of time are indicated as 
solid/dashed lines in the insertion of fig. 6.4(b) for T = 150, 200 and 300°C and as 
dashed/solid lines in the inset of fig. 5.3(a) for T = 250 °C. 


T(°C) 

150 

200 

250 

300 

P 

0.51(4) 

0.41(5) 

0.24(4) 

0.21(5) 

1 /z 

« 0.07 

0.6(1) 

0.62(2) 

0.7(1) 


Table 6.3: Values for the universal exponents f3(T) and 1 /z(T) for surface fluctuations 
of CdTe grown on Si(001) substrates by HWE with F « 2.2 A/s within the range of 15 
to 240 min and for several deposition temperatures. 


6.4 Universal Distributions 

In this section, we supplement our previous studies with a deep analysis of height, squared 
roughness local and maximal relative height distribuitons. As we shall demonstrate, this analy¬ 
sis, rather than a complementary one is, in some cases, essential for unveiling the Universality 
Class of a given growth. 
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6.4.1 CdTe surface fluctuations at T = 150 °C: Poissonian Growth 

Figure 6.7 shows the rescaled height distributions for different growth times for CdTe surfaces 
grown at T = 150 °C. The HDs arc compared to the Gaussian distribution, expected for un¬ 
correlated growing surfaces. One notices that experimental data collapse very well with the 
HD-Gaussian near of the peak (inset (A) of the fig. 6.7) as well as at the tails - in at least four 
decades around the peak. 



[h - <h>]/o h 


Figure 6.7: Rescaled height distributions at mean null and unitary variance for dif¬ 
ferent growth times (symbols) from CdTe surfaces grown at T = 150 °C. The experi¬ 
mental data are compared to the Gaussian (solid cyan line). Inset (A) highlights the 
very good collapse between experimental data and the Gaussian close to the peak in 
a linear x log plot. Insertion (B) exhibits the experimental S and K evolving in time. 
Solid/dashed cyan and magent lines refers to the S/K values for the Gaussian-HD and 
KPZ-HD, respectively. 


Experimental values for the Skewness (S) and Kurtosis ( K ) are displayed in the insertion 
(B) of the fig. 6.7, where their values seem to fluctuate around zero, as expected from interfaces 
having a Gaussian-HD. Adding to this result those coming from the scaling roughness analysis 
(see table 6.3) as well as those from the local dynamic, one can conclude that, at this low 
temperature, height-field fluctuations of CdTe films are Poissonian [22,159]. 

The S and K values as function of the box size (I) for SLRDs and MRHDs also go 
in favor to a uncorrelated growth - see tigs. 6.8(a) and (b). As discussed in the sections 3.2 
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Figure 6.8: S and K values for (a) SLRDs and (b) MRHDs from CdTe surfaces grown 
at T = 150°C by 15 min (orange left triangles), 30 min (blue up triangles), 60 min 
(green diamonds), 120 min (red squares) and 240 min (black circles). No one con¬ 
vergence for S(l) and K(l) curves are observed, as expected for uncorrelated growing 
surfaces. 


and 3.3, these “box” distributions must be evaluated in the range of £ <C l <C £. When l 
satisfies this limit, S and K values become independent of l and converge to universal values 
[129, 130,134, 136]. However, when the surface is uncorrelated, £ keeps very small being of the 
same order of £ (when there are grains at the surface). Thus, no one convergence of S and K 
should be observed, as can be confirmed in the figures 6.8(a) and 6.8(b). 

6.4.2 Random-to-KPZ Crossover in CdTe Surface Fluctuations at T = 200 °C 

HDs for CdTe films grown at T = 200 °C can be seen in the fig. 6.9. Very interesting, one 
notices that the distributions at initial growth times are very far from the KPZ one. Instead, 
they are closer to a Gaussian. This is confirmed by S(t) and K(t), which are closer to zero 
than to the KPZ values - see inset B of the fig. 6.9. At long times, however, S and K values 
seem to converge to the KPZ ones and a very good agreement with the universal KPZ-HD is 
observed. This disagreement/agreement with a KPZ-HD is reinforced taking into account the 
collapse near of the peak (see inset A of the fig. 6.9) for surfaces grown at t = 15 min and 
t = 240 min, respectively. 

Despite of large error bars related to S and K values, which stem from low statistics, the 
agreement between the longest growth time available and the expected KPZ values (see table 
3.1) is remarkable, namely S /-240 = 0.43(5) and K i= 24 o = 0.5(2), whereas KPZ values are 
S = 0.42(1) and K = 0.34(2). Based on our previous local and roughening scaling study, 
we can argue that the HDs above show a crossover in time toward the KPZ regime. In 
particular, this kind of crossover is widely studied numerically (see [163-165] and ref. therein) 
and this experiment, as far as we know, is the first and clean experimental evidence of such 














6. Temperature Effect on CdTe Growth Dynamic 


58 



Gaussian 

— 

KPZ 

◄ 

15 min 

▲ 

30 min 

♦ 

60 min 

■ 

120 min 

• 

240 min 


240 min 


(B) 


' . 1 

■ 

■ 

H, - 

... 



** 

: £ 

rm 


0 60 120 180 240: 

t (min) 



-2 0 2 

[h - <h>]/o h 


Figure 6.9: Rescaled height distributions, at mean null and unitary variance, for dif¬ 
ferent growth times (symbols) from surfaces grown at T = 200°C. The experimental 
data are compared to the Gaussian (solid cyan line) and to the numerical KPZ curve 
(solid magenta line). Inset (A) highlights the very good collapse between experimental 
and numerical data close to the peak in a linear x log plot. Insertion (B) exhibits the 
experimental S and K evolving in time. Solid/dashed cyan lines refers to the Gaussian 
S/K values, whereas the magenta line indicates the KPZ S/K values. 


random-to-KPZ growth in d = 2 + 1. 

The SRLDs and MRHDs also support these findings, as exhibited in the figures 6.10 and 
6.11. At initial times, SLRDs curves do not present the stretched exponential decay at the right 
tail, which is known to be a KPZ landmark [132]. Moreover, S and K values for these times do 
not converge to the KPZ one as the box size is increased (see the fig. 6.10(b)). The comparison 
between the SLRD for surfaces grown at t = 15 min and the KPZ-SLRD around the peak gives 
more evidences of such discrepancy - inset of fig. 6.10(a). Nevertheless, the situation is totally 
different at long times, where the agreement between the experimental data and the theoretical 
one reaches more than four decades around the peak with S and K values converging to the 
KPZ values as the continuous limit is approached, constrained to ( <C l <C £. Results coming 
from the MRHDs (fig. 6.11) are not so accurate as those emerging from SRLDs, but they also 
corroborate that, for CdTe thin-films grown at T = 200 °C, a Random-to-KPZ crossover takes 
place. 
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Figure 6.10: (a) Rescaled squared local roughness distributions for samples grown at 
T = 200 °C, by 15 min (orange left triangles), 30 min (blue up triangles), 60 min (green 
diamonds), 120 min (red squares) and 240 min (black circles). The curves shown are 
those whose the box size is indicated below the legend. They are different for each 
t because the appropriate interval C < l < £ depends on t. The inset shows the 
poor/nice collapse around the peak for the thinnest/thickest film grown, (b) Skewness 
and Kurtosis values as function of the box size l. 
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Figure 6.11: (a) Rescaled maximal relative height distributions for samples grown at 
T = 200 °C, by 15 min (orange left triangles), 30 min (blue up triangles), 60 min (green 
diamonds), 120 min (red squares) and 240 min (black circles). The curves shown are 
those whose the box size is indicated next to the legend. They are different for each t 
because the appropriate interval C < ( < £ depends on t. The inset shows the good 
collapse around the peak for the thinnest and the thickest film grown, (b) Skewness 
and Kurtosis values as function of the box size l. 
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6.4.3 KPZ Growth with Deposition Refused: CdTe Surface Fluctuations at T = 

300 °C 

Finally, the analysis of films grown at T = 300 °C enclose this section and bring to out the real 
UC of CdTe thin films grown at high temperatures. Figure 6.12 depicts the experimental HDs 
being compared to numerical KPZ, MFI (linear-MBE equation) and VLDS (non-linear MBE 
equation) curves. 



[h - <h>]/o h 

Figure 6.12: Rescaled height distributions at mean null and unitary variance for dif¬ 
ferent growth times (symbols) from surfaces grown at T = 300°C. The experimental 
data are compared to the numerical KPZ curve (solid magenta line) considering A < 0, 
to the universal HD-MH (solid cyan line) and to the universal HD-VLDS (dashed green 
line). Inset (A) highlights the very good collapse between experimental and KPZ nu¬ 
merical data close to the peak in a linear x log plot. Insertion (B) exhibits the ex¬ 
perimental S and K evolving in time. Solid/dashed lines refers to the expected S/K 
values. 


Just to mention, we have checked the universality for the HDs of these diffusion-dominated 
classes by integrating the linear equation as well as simulating several discrete models (see 
[94,96,166] and ref. therein) that belong to the MH and VLDS classes. As an additional note, 
the KPZ HD plotted is that one for A < 0, which is the same HD plotted in the previous figures, 
but reflected around the origin. Anyway, one can observe that HDs for early growth times do 
not agree with distributions of any class, but, for long times, once again, a reasonable agreement 
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with the KPZ-HD is observed. The left tails seem to converge faster than the right ones to the 
KPZ curve - see HDs for t > 30 min. As demonstrated by the AFM images (fig. 6.1(e) and (f)) 
and by the local exponent a \ (inset of the tig. 6.4(c)), the top of grains become smoother as the 
time evolves. It leads the right tail of the HDs to move in the left direction and, consequently, to 
approach to the KPZ-HD. 

Unlike the exponents, the experimental HDs discard the MH and VLDS classes as pos¬ 
sibles UC of the CdTe growth. For the largest time available, the inset A of the fig. 6.12 presents 
a nice collapse around the peak between the experimental result and the numerical one. More¬ 
over, it is clear seen that S(t) and K(t) values are approaching to the KPZ ones, within the error 
bars - inset B of the fig. 6.12. Particularly, for t = 240 min, one has Si= 24 o = —0.2(2) and 
Kt =240 = 0.3(2), which agree with the KPZ values (table 3.1), for A < 0, considering the large 
eiTor bars. 



Figure 6.13: (a) Rescaled squared local roughness distributions for samples grown at 
T = 300 °C by 15 min (orange left triangles), 30 min (blue up triangles), 60 min (green 
diamonds), 120 min (red squares) and 240 min (black circles). The curves shown are 
those whose the box size is indicated below the legend. They are different for each t 
because the appropriate interval C < l < £ depends on t. The inset shows the good 
collapse around the peak for the thinnest and the thickest film grown, (b) Skewness 
and Kurtosis values as function of the box size l. 

To gain further evidence of KPZ scaling, the SLRDs and MRHDs have been calculated 
and are depicted in the figures 6.13 and 6.14, respectively. From the left tail and from the peak of 
the SLRDs (fig. 6.13 and its inset), one notices that experimental results are very well described 
by the KPZ curve, whereas the decay at the right tail can not make a clear distinction between 
the KPZ and the MH class. In the fig. 6.13(b), S values converge to the KPZ ones, in the 
meantime that I\ values also can not be distinguished from that expected for the KPZ-MRHD 
and MH-MRHD. In this situation, where the exponent analysis point to KPZ or to MH growth, 
the MRHDs are not suitable to make a separation between these two classes. As showed in the 
fig. 6.14, the expected S and K values for the KPZ and the MH class are very close between 
each other and, from the experimental point of view, one can not distill the UC based only on 
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Figure 6.14: S and K as function of the box size for MRHDs for samples grown at 
T = 300 °C. Solid/dashed line indicates the expected values for the S and K for the 
KPZ (magenta line) and for the MH class (cyan line). 


this result. 

These findings help us to see how hard is determining the UC of the two-dimensional 
experimental films. Anyhow, guided by the local-structure evolution, roughening scaling results 
& distributions analysis, we can conclude that the CdTe interface for T = 300 °C evolves 
according the KPZ equation (eq. 3.1) with A < 0. This indicate the existence of a mechanism 
of refusing the deposition of particles depending on the local slopes as occurs in the famous 
RSOS model [128]. A possible explanation for A < 0 is that the sticking coefficient is smaller 
in regions with very large slopes at surface. Indeed, one notices in the fig. 6.4(b) that ((Vfi) 2 ) 
for T = 300 °C is very larger than that for lower T. This can explain why this effect appears 
only at high T [61,167]. 


6.5 Final Remarks 

A detailed study on how the deposition temperature affects the CdTe growth dynamic in a wide 
range of temperature [150, 300] °C has been performed. The local mound evolution has been 
explained in terms of an interplay between a local energy bander Eqb that tends repelling the 
diffusion and the deposition of particles at defect sites (yielded by the colision of GBs of neigh¬ 
boring grains) and the effectiveness of the relaxation process at these GBs, which tends to elim¬ 
inate defects on the active zone of the CdTe interface. In turn, the last mechanism depends 
sensintively on the substrate temperature. Thus, for T = 150 °C, the low diffusion at the GBs 
of colided neighboring grains prevents coalescences and the propagation of correlations at in¬ 
terface, so that inter-grain fluctuations are described by Possianian process. This dynamic also 
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happens at short times for films grown at T = 200 °C, however the relaxation at (and around) 
the GB's of colided neighboring grains gives rise to coalescences and, hence, an asymptotic 
correlated growth. For higher T, where diffusion rate is higher, these processes start early as 
well as the KPZ scaling. 

In terms of KPZ equation (eq. 3.1), the random growth at low T implies v « 0 and 
A « 0. For T = 200 °C, one expects A > 0, but small, so that growth is dominated by 
noise/non-linear effects initially/asymptotically. The absence of a crossover when T = 250 °C 
[63] indicates that A is positively large. The same inteiplay also takes place in this temperature, 
but as the diffusion is still more operative than that for lower T, coalecence processes start 
early enhancing the activity the grain-packing (discussed in the previous chapter) and, hence, 
the strength of A > 0. Thus, A(T) seems to be a positive increasing function in this range of 
T. However, for T = 300 °C, the evidences of KPZ scaling and the negative skewed HDs- 
KPZ reveals A < 0, which is typical of KPZ systems where there is deposition refuse as in the 
RSOS model [128]. A possible explanation is that the the sticking coefficient becomes smaller 
in regions with very large slopes at surface. This reasoning is corroborated by ((V/i) 2 ) for films 
grown T = 300 °C, where one confirms the presence of larger slopes than those for lower T. 
This can explain why this effect appears only at high T [61,167]. Anyway, it is astonishing that 
so contrasting KPZ mechanisms can emerge in the CdTe growth system only adjusting T. 




Conclusions and Perspectives 


In this work one has performed a detailed study on the growth dynamic at both short- and large- 
wavelength scales of CdTe thin films grown on Si(001) substrates by Hot Wall Technique. The 
results have provided, for the first time, a clear and robust evidence of a system in d = 2 + 1 
dimensions which belongs to the Kardar-Parisi-Zhang class. Here, the phrase “robust evidence” 
means that the prove goes beyond the comparison with critical exponents, as confirmed by 
the (rescaled) height distributions, squared local roughness distributions and maximal relative 
height distributions. In the meantime, this work demonstrates the universality of these KPZ 
distributions, giving them a reliable reality beyond numerical simulations. 

Along the analysis, one has found several pitfalls hampering the extraction of asymptotic 
scaling exponents. For instance, through the power-law in the wi oc x 1 plot, it was not possible to 
unearth the roughness exponent. Rather than, only the geometrical exponent, ai, which does not 
provide any information about the Universality Class (UC) of the system, was found. Another 
difficulty is the relation between the average mound size (Q and the correlation length (£). The 
validity of C ~ should be taken as true only for an appropriate range of time, for which local 
peaks on the top of large mounds do not strongly underestimate the measurement of Thus, 
an initial deep-visual inspection of surfaces imaged by some microscope technique becomes 
important. 

The KPZ mechanism occurring in CdTe films grown at T = 250 °C is related to the 
form in which new particles attaches to grain boundaries (GBs) of colided neighboring grains, 
do not filling all available space in their neighborhood and, thus, generating excess of velocity, 
the KPZ landmark. As will be discussed below in details, this packing mechanism also should 
be present in other temperatures. Due to inexorable experimental obstacles it was not possible 
to validate the KPZ Ansdtze in d = 2 + 1 dimensions. However, this is a natural extension of this 
work. 

The effect of the deposition temperature (T) on the growth dynamic has also been dis¬ 
tilled in broad range of T, namely, T £ [150, 300] °C. A relation between short- and large- 
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wavelength dynamics has been established. One has been found that the mound evolution is 
dictated by the interplay between the formation of defects at grain boundaries of colided neigh¬ 
boring grains and the relaxation process induced by diffusion and deposition of particles toward 
these regions. A simple Monte Carlo model corroborate this reasoning. This interplay leads 
to different scenarios at large-wavelength fluctuations as T increases. For T = 150 °C, the 
low diffusion at the GBs prevents coalescences and the propagation of correlations at interface, 
so that inter-grain fluctuations are described by a Poissonian process. Scaling exponents and 
distributions give support to these reasonings. 

For T = 200 °C, however, a more complex scenario has been found. In the range of 
growth time studied, the relaxation process overcomes the effects of the barrier at (and around 
of) the GBs and gives rise to a crossover in the CdTe growth dynamic. Results from dynamical 
scaling and distributions also corroborate with the presence of this crossover. Roughness scaling 
analysis, nevertheless, have not been able to, convincingly, point the UC of the CdTe growth at 
this deposition temperature. The supplemental study based on the distributions, in turn, have 
proved itself be essential for unveiling the UC of this system, once weaker finite-time effects 
than that in the standard roughness scaling has been found. Thus, it was possible to find a 
Random-to-KPZ crossover taking place in the CdTe growth, at the same time that the first robust 
experimental realization of such crossover in two-dimensional systems is demonstrated. 

The growth of CdTe films at T = 300 °C carries also its distinct importance. Although 
dynamical scaling analysis do not allow making a clear distinction among diffusion-dominated 
equations and the KPZ one for describing CdTe fluctuations at this high deposition temperature, 
the distributions discard the formers and point out the presence of KPZ growth with A < 0, 
differently of lower temperatures. This unanticipated result has been related to the decreasing 
of the sticking coefficient in regions with very large slopes at surface. The reasoning is cor¬ 
roborated by the squared-local slopes, which arc about five times larger than those for lower 
temperatures, as can be noticed in the Atomic Force Microscope images. This can explain why 
this KPZ mechanism appeal's solely at high temperatures. 

In terms of the KPZ equation (eq. 3.1), the results reveal that, for T = 150 °C, the 
“surface tension” (y) and the excess of velocity (A) are very near of zero, so that noise dominates 
the growth. For CdTe films grown at T = 200 °C, however, the Random-to-KPZ crossover 
found indicates A > 0, but small, so that non-linear effects overcome the noise only at long 
growth times. In turn, a clean evidence of the KPZ scaling since inital growth times is present 
in CdTe films grown at T = 250°C, in which A is positive and large. At the moment that all 
conspire to conjecture A(T) as an increasing function of T, the results coming from CdTe films 
grown at T = 300 °C show that A is negative. This suggests that i) is possible to adjust T at some 
specific Tew > with T <G ]250, 300[°C, in order to obtain A = 0, i.e, to have a growth described 
by the Edwards-Wilkinson equation and ii) to set T^ew-kpz) = Tew ± ST, so that EW-to- 
KPZ crossovers emerge into the dynamic of CdTe surface fluctuations. In summary, it is possible 
to adjust the KPZ non-linearity in the CdTe system only adjusting the deposition temperature. 
Figure 7.1 summarizes theses discussions, making clear different regimes that takes place in 
the growth as T is increased, namely: region I refers to the temperature interval in which a 
Poissonian(Random)-to-KPZ crossover rules the dynamic. This first region is localized at the 
rigth border of temperatures close to 150 °C, where A —> 0 and where the Poissonian(Random) 
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Figure 7.1 : Conjecture for the behavior of the excess of velocity, A, as function of 
the deposition temperature T in the CdTe/Si(001) system. Circles refer to experimen¬ 
tal points, and dashed lines indicates expected different regimes (I to IV) ruling the 
coarsening-grained dynamic. The “x” point denotates temperatures for which global 
fluctuations are expected to belong to Random (green x) and EW (blue x) Classes. 


growth emerges; in the region II a clean KPZ scaling emerges and precedes the thin III regime, 
limited by T = Tew ± ST and symmetric in relation to the Tew point, in which the pure 
EW growth (ST = 0 °C) and the EW-to-KPZ crossover are expected to take place. Finally, at 
high temperatures (although expected do not be much higher than 300 °C), the KPZ scaling is 
recoreved with A < 0. 

We finish this work truly believing that the noveil procedure for investigating the UC of 
growing surfaces open and will motivate (as exemplified by T. Halpin-Healy and G. Palasantzas 
[64]) a new confident perspective in the field, as well as the application of these methods in 
previously studied systems. Our system also offer the experimental possibility of a close study 
of Random-to-KPZ and EW-to-KPZ crossovers. 












More Details About Continuum Growth Equations and 
Universality Classes 


A.l Random Growth equation 

Concerning on EW and linear-MBE equations, one can easily see the particular case, where v 
and K are nulls respectively, resulting in the so called Random Growth equation: 


d t h = rj(x,t). (A.l) 

Once just the noise rules the growth dynamics, there arc no correlations to be propagated 
through the system and we can argue z —> oo. The a exponent is also ill defined, since the 
system never reaches the stationary regime. Nevertheless, we still can find as the variance of 
such random interfaces evolves with the deposition time. The result is anticipated if one looks 
at these surfaces as being composed by a set of independent Brownian motions. 

Integrating the random eq. one obtains: f dth(x,t)dt = f r/(x, t)dt. Averaging both 
sides, one has ( h(x,t )) = 0. But, if instead we square and average, one reaches (h 2 (x,t.)} = 
2 Dt. So, the variance grows as in the eq. A.2, defining the “growth exponent” 1 /3 = 1/2. 

w 2 {t ) = ( h 2 ) - ( h) 2 = 2 Dt. (A.2) 

1 In accordance with our previous definition, j3 = a/z. However, the random growth is a particular case, where 
one can find a power-law in time for the variance, even when the spatial- and temporal-statistical-fractality exponents 
are ill defined. In this sense, the j3 exponent is fundamentally different front those previously cited. 


68 



Appendix A - Details of Continuum Growth Equations and Universality Classes 


69 




Figure A.1: Schematic of the chemical potential depending on the local curvature. 
At left to rigth: positive curvature, showing a valley, where there is a local minimum 
of /x(x,t) because particles sticked at the bottom of the valley will not move away 
from there. On a terrace R oo and n{x,t) is null everywhere. Negative curvature, 
showing an island, where n(x,t) has a local maximum value because particles on the 
top of the island break off their bonds and hop to a closest neighbor site where the 
free energy of the interface is decreased. 


A.2 A derivation for the Linear-MBE equation 

It is well known that in MBE environments, under usual conditions, desorption 2 is a negligible 
process and, hence, the therm V 2 /r in the eq. 2.19 should not play an important role. Formation 
of overhangs and bulk defects are uncommon in the ideal epitaxy [155, 156] and discard the 
non-linear (V/i) 2 term. Our starting point for building a continuum growth equation for MBE 
environments is the conservation law: 


d t h = -V • J + r)(x,t), (A.3) 

where J represents a parallel current of particles diffusing onto the interface. 

As interfaces growing in MBE chambers are driven by chemical bond in order to min¬ 
imize their surface free energy through diffusion, J should be oc —V/i. But the simplest as¬ 
sumption, in turn, points out that the chemical potential is proportional to — 1 /R, where R is the 
local curvature - see fig. A.l. The bottom of a valley (positive curvature, R > 0), for instance, 
is the site where a particle have major number of neighbors and move away from there will be 
difficult. So, the bottom of a valley is a minimum for The opposite happens for the 

topest site of a island (R < 0), where there is a local maximum for fi(x,t). On a terrace, all 
particles have the same coordination number, and locally fi = 0 everywhere. Hence, on these 
considerations one can argue straightforwardly /i(x. t) oc — V 2 /i. Consequently, J oc V(V 2 /i) 
and one reaches to the eq. 2.23, after eq. A.3. 

2 See Kinetic Effects in the Appendix section C.4 
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A.3 The Non-Linear MBE equation and the VLDS class 

In models with conserved dynamic, there is necessity of inserting non-linear terms, since a > 1 
breaks the hypothesis that continuum approximations are valid in the limit of small slopes, 
Vh <C L“ -1 [22], As pointed out by Wolf and Villain [72], the lowest non-linear term obeying 
conserved dynamic is the therm V 2 (V/i) 2 . The geometrical interpretation (see ref. [168]) indi¬ 
cates that particles landing at high steps (large derivatives) relax to lower steps (small deriva¬ 
tives), corresponding to the “high-temperature” regime in MBE environments, where atoms 
break their bond and diffuse, with larger bias, to small slope kinks. The non-linear MBE equa¬ 
tion (eq. A.4) was proposed formally in 1991 by Lai and Das Sarma [168]. 


d t h(x, t ) = -/TV 4 /?, + AiV 2 (V/r) 2 + r/(x, t), (A.4) 

Ai has the same dimension of v and accounts for the strength of diffusion of particles at local 
slopes towards small derivatives. 

Critical exponents related to the non-linear MBE equation can be found by using Renormalization- 
Group (RG) approaches [22,168]. The method predicts: 


a = 


4 - d s 
3 



(A.5) 


Other very important result coming from RG analysis is the hyper-scaling relation valid 
for any growth with conserved dynamic and non-conversed noise'. 


z — 2a — d s = 0 


(A.6) 


The set these exponents compose the VUlain-Lai-Das-Sarma (VLDS) Universality Class. 
There are a small number of discrete models belonging to the VLDS class, among which pos¬ 
sibly are the conserved restricted-solid-on-solid (CRSOS) model, proposed by Kim et al. [160] 
and the Das Sarma and Tamborenea one (DT) [91]. A very good numerical study and descrip¬ 
tions of these models can be found in the Ref. [166]. 






Anomalous Scaling 


Anomalous scaling occurs whenever a surface is not self-affine [67]. It means that local and 
global fluctuations do not evolve in the same way. Initially, the therm anomalous was used to 
describe the growth of interfaces in which the exponent a was found to be larger (or equal) 
than 1 [66, 169-171], Indeed, this comes from the fact that, if a ^ 1, there is a divergence 
of fluctuations in the steady state, making the surface super-rough, as indicated by the Family- 
Vicsek ansatze (eq. 2.11): W sa t/L ~ L a ~ * l . 

The current situation, however, is very different from that of almost twenty years ago. 
Nowadays, one knows that several systems display anomalous scaling, independently of the a 
value (see [70,71] and references therein). Moreover, there are more than one possible scaling 
form to describe the local-fluctuations dynamic of a growing interface, being the FV one only 
a particular case where there is no difference between the local and global one [70]. A wide 
number of experimental studies have also corroborated this conjecture, among which are the 
development of wood fractures [172], the electrodeposition of Cu [158], the dissolution of pure 
iron [85] and the deposit of colloidal particles at the edge of evaporated drops [173]. 

The initial clues of a break in the FV hypothesis were found in the one-dimensional 
Wolf-Villain model [66, 73] and, subsequently, in other linear and non-linear growth models 
dominated diffusion [73,96, 170, 171], From there, analytical and numerical treatments were 
strongly devoted to this issue, leading to the conclusions: 

i) The condition a < 1 is not sufficient to prevent the system of exhibiting the anomalous 
scaling [67]. 

ii) There exist three types of anomaly: super-rough, intrinsic [68] and faceted [70] (they are 
defined below). 

iii) Non-conserved growth models do not exhibit anomalous scaling and, in particular, intrin¬ 
sic scaling can not occur in local growth models [69,71]. 
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It is worth mention that the presence of the anomaly affects only local correlation func¬ 
tions, keeping the global roughness, for instance, scaling according to the standard FV scal¬ 
ing [70] - this is proved in the following. For now, consider the Fourier transform of the function 

h(x, t ) as, 


h(k, t) = L ds / 2 ^[/t(x, t) — (/i(i))] exp(zk • x), (B. 1) 

X 

where k arc the wave-numbers surviving within the space bounded by the substrate. 

The called structure factor [S(k, t) ] or power spectrum (eq. B.2) is a correlation function 
measuring the fluctuations in the reciprocal space. 


S(k, t) = (h(k, t) h(-k, t)) = (|h(k, t )| 2 ). (B.2) 

S(k. t) can be related to the global roughness and to the height-height correlation func¬ 
tion ( Ch ) by the eqs. B.3 and B.4, respectively, 

1 _ f rj^ s 1c 

W\L,t) = W- Y,S(k,t) = J ^S(k,t), (B.3) 

C h (l,t) k /[l-coa(k -x)] p * (B.4) 

where the integrals run within the interval 2 tt/L < k < 27r/a for each k direction, with a being 
the lattice parameter. 

The basic tools required to derive the different forms that anomalous scaling can appear 
are written in order in the equations above. Inserting the FV ansatze (eq. 2.11) into the relation 
of eq. B.3, it is straightforward showing that the power spectrum behaves as: 


S(k,t) = k-( 2a+d ^s FV (kt 1/z ), (B.5) 

where: 


S F v(u) ~ 


const, if u» 1, 
u 2a+d , if u <C 1. 


(B-6) 


It means that for short-length scales (k S> 1/0, the power-spectrum is time-independent 
and the curves for all times show a power-law scaling with k-(' 2fi + d ») . xhis is the usual FV 
scaling in the reciprocal space. But now, look what happens when this scaling relation is inserted 
into the eq. B.4. When one sets a > 1 , the integral in eq. B.4 becomes divergent in the limit of 
l <C £ for L —» oo, and a —> 0 [68]. Taking the limit l -C £ first, and keeping L and a fixed, one 
obtains a different scaling relation for Ch- 


C h (l,t)~ 


Z 2 t2(«-1 )A if l < £ < L, 
l 2 L 2{a - l \ if ^<L«£. 


(B.7) 


So now a different roughness exponent has emerged and, locally, Ch(l,t ) /2(m oc ) ? 

where cq oc f a and cq oc = 1. Moreover, a dependence in time for Ch at wavelengths l <C £ 
appears in the system scaling with t 2K , where: 
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k = (a- Otloc)/ z. 


(B.8) 


Translating this scaling for the context of the fig. 2.3, one would see the curves for different 
growth times, shown in that figure, shifted to up as the time evolves for l <C £. Figure B.l(a) 
shows very clearly this non-usual behavior by calculating the square-local-slope evolution (eq. 
B.9) for the one-dimensional Wolf-Villain model [66]. 


((Vh) 2 )~C h (l = l,t)~t 2K . 


(B.9) 





Figure B.1 : Example of anomalous scaling. In (a) one has the function G(1 ,t) (corre¬ 
sponding to the C h (l,t) function) for several systems sizes L (from the bottom o the 
top L = 10, 15, 30, 60, 120, 240, 800 and 1024, respectively), (b) Shows the non¬ 
usual behavior for the structure factor S(k,t). From the bottom to the top the curves 
refer to growth times t = 2 4 , 2 7 , 2 10 , 2 13 , 2 16 , 2 19 , 2 22 , respectively. Figure extracted 
from ref. [66]. 

In the WV model, it was found k = 0.19(1) during a transient regime, limited by both 
finite-size and finite-time collections [66]. But this feature of non-usual behavior of correlations 
functions at local scales is not the complete scenarium in which the anomaly can emerge. Indeed, 
models incorporating random diffusion [68] and diffusion-dominated [174] presented a non¬ 
usual behavior not only for Ch(l, t ), but also for the power spectrum, i.e [S (k l/£,t)] is time 
dependent - see fig. B.l(b). Thus, it became clear that the FV ansatze written in the reciprocal 
space (eq. B.6) could not be the general scaling dictating local fluctuations in the surface growth 
context. 

The generic dynamic scaling in kinetic roughening has been developed along of many 
years [67-69], but a formal deduction has been achieved in the year of 2000, where Ramasco et 
al. [70] have proposed a scaling for s(u ) as, 



u 1, 
«< 1 , 


(B.10) 
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where a s is a new exponent called spectral roughness exponent. 
Thus, the power-spectrum, genetically, scales as: 


s(M) ~ 


fc -(2a.+d.) j if k^l/t 1 ^, 
t * , */ k « 1/T 1 / 2 . 


(B. 11) 


Inserting eq. B.10 in eq. B.3, one proves that global fluctuations are dictated by the 
standard Family-Vicsek scaling (eq. 2.11), independently of a and a s values. Nevertheless, 
as we have shown above, the local scaling is modified, depending on the order of the several 
limits involved and on values of the exponents [70]. In a general way, two major cases can be 
distinguished, namely, a s < 1 and a s > 1. For the former, one obtains the same that in the eq. 
2.10, but with L —> l and a —> a s , i.e: 


W(l,t) = f as< \(l /£) , 


(B.12) 


with, 


fa 3 < lO) ~ 


u“ s , if u <C 1, 
const, if m»1, 


(B.13) 


implying that a s = cq oc [67,68]. This is called intrinsic anomaly: one has the power-spectrum 
and w(l, t ) [ Cjfl, t) | behaving non-trivially. Family-Vicsek scaling is recovered when a = cq oc 
[70], 

Otherwise, for a s > 1, the integral in the eq. B.4 is divergent when L oc. Keeping L 
fixed, one finds that: 


fa s >l(u) ~ 


u, if u <C 1, 

const, if u 1. 


(B.14) 


This implies that cq oc = 1, independently of a s . Additionally, whether a = a s , so the 
power-spectrum evolves trivially in time, whereas w(l, t) and Ch(l, t) do not. This is the super- 
roughening anomaly, discussed at the beginning of this section. Flowever, when the equality 
between a and a s is not fulfilled, so a new anomaly stems, namely, the faceted scaling [70]. 
The main characteristic of the faceted anomaly is that it can detected only by using the power 
spectrum because there is not a constrain between a and cq oc , i.e, these exponents can be equal 
(k = 0) or different (k / 0). 

The summary of all conditions leading to the different scaling which can appear in sur¬ 
face growth is, hence: 


7/ ot s < 1 ^ Q7oc — 

< 

7/ ot s > 1 cq oc — 1 


J a s = a =y> Family — Vicsek , 
| a s / a =y> Intrinsic, 

( a s = a => Super — rough, 

\ a s yl a => Faceted. 


(B.15) 


From the experimental point of view, several studies have found the different anomalous 
scalings. For instance, a transition from intrinsic to faceted scaling has been reported during 
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dissolution of iron [85], whereas a faceted dynamic has been claimed in the growth of CdTe 
on glass substrates [15]. Intrinsic anomaly seems to appear in the growth of TiN thin films by 
reactive sputtering [175] and during the deposit of colloidal particles at edge of evaporated water 
drops [173]. Other studies involving anomaly can be found in the Ref. [65]. 

Progress has also been made in the sense of classify when a possible anomaly can appear 
depending on the continuum equation considered. Lopez et al. [71] has used a transformation 
T = V/i to study the roughness (VFx = ((V/i) 2 )) of surfaces described by slopes fields. When 
such surfaces are rough, it implies k > 0 and one has anomaly (disconsidering the faceted case). 
From this starting point, one can show that non-conserved growth models, such those ruled by 
the Kardar-Parisi-Zhang (KPZ) class 1 can not exhibit anomaly, whereas local growth models or 
display FV scaling or super-roughening. They can not be exhibit intrinsic anomaly [71]. 


'See chapter 3. 





Concepts for Surfaces in Equieibrium and for the Thin 
Film Growth 


C.l Surface Tension and Equilibrium Shape 

A system is in equilibrium when its macroscopic properties do not change appreciably with the 
length scale of measurement time [80]. In this regime, there is a suitable set of thermodynamic 
coordinates describing the system, which includes generalized displacements {x} and their re¬ 
lated generalized forces {Q} [80,81]. Examples for the former are the volume (V) for a gas or 
the area {Ay) for a film, whereas their related forces are the pressure (-P) and the surface tension 
( 7 ), respectively. 

Depending on what coordinates are fixed in an experiment, the system state can be 
specified by appropriate thermodynamic potentials. For instance, for processes approaching 
to the equilibrium isothermally [the temperature (T) is constant] without work (d W = 0), the 
Helmholtz free energy Fh = Fh(T, {x}) is well defined. Otherwise, if only chemical work 
(&W c he = Yli HidNi, where /i, is the ith-chemical potential and N t is the particle of ith-nature) 
is absent, one can use the Gibbs free energy G = G(T, {Q}). Both are minimum at the equilib¬ 
rium and contains all thermodynamic information accessible to the system [80]. 

Once one dealing with surfaces (or interfaces), one shall consider that boundary atoms 
are more energetic than those in the bulk because they have dangling bonds. Thus, for creating 
a new surface of area dAf, one must provide an amount of energy (IFh\t,v = ydAj neces¬ 
sary for breaking chemical bonds. In this situation, 7 is also interpreted as the excess of free 
energy per unit of area, often called superficial energy [2]. Otherwise, if there arc particle migra¬ 
tions between the bulk and the surface, appears a term of chemical work and one must provide 
dr\T,v = 7 dAf — Y, t N,df.i t , with S being the grand potential. 

Anyway, the important is: a surface minimizes its energy decreasing its area or changing 
7 . The first case is the reason for which a water drop does not spread on the dew, and for 
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which without gravity its shape is perfectly spherical 1 . Cappilarity phenomena and soap bubble 
dynamics are also dictated by the surface tension [2,3,176]. 

Typically, metals have higher 7 than oxides 2 and organic structures [1]. Water liquid- 
vapor interfaces has 7 = 72.94 niN/m, whereas that ones formed by glycerine or methanol have 
48.09 and 22.50 mN/m at 20 °C [176]. If the surface cannot decrease its area, 7 variations leads 
to modifications in the atomic arrangements and arc responsible for surface reconstructions in 
crystals [1-3]. 



Figure C.1: SEM image of Pb crystal at equilibrium viewed from the [011 ] azimuth for 
(a) 300 °C, (b) 320 °C and (c) 327°C, above the Pb melting point. Extracted from [2], 


The shape of an object at the equilibrium can be predicted by variational principle of the 
suitable thermodynamic potential [2]. For a crystal, the surface tension depends on what face is 
considered and 7 = 7 (hkl, T ), where hid are the Miller’s indices. This anisotropy governs the 
crystal shape and, as the temperature increases, it tends to be minimized [3]. Above the melting 
point, the anisotropy vanishes and the material must recover the spherical form [2,3]. In the fig. 
C. 1 one can seen equilibrium shapes for Pb crystals at different temperatures. As temperature 
increases the anisotropy decreases and fewer facets are formed. Above the melting point, where 
Pb is liquid, the spherical shape is achieved. 


C.2 Nucleation 

Nucleation is a fundamental process happening during a phase transition [1], It is important 
for thin films because the structures appearing in the submonolayer regime strongly affects the 
dynamic of the growth [22,159]. Simple models for nucleation are based on the liquid drop 
model, which also works for solids. The basic difference is that solids have two additional 
features: elasticity and commensurability [3]. These effects will be discussed along this section. 
Nucleation models are important just as a qualitative fashion once one deals with few atoms 
and, hence, a continuum description (on thermodynamic basis) of the phenomenon becomes 
doubtful [ 1 ], 

'The spherical shape is that one that minimizes an area for a given volume, since 7 is isotropic [176]. 

2 This is intuitive, once metallic surface tends to become oxided. 
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Deposition 



Figure C.2: Schematic of nucleation processes in a gas-to-solid phase transition. 


Consider a possible heterogeneous nucleation of a condensed film on a crystalline sub¬ 
strate (see a scheme in the fig. C.2). Deposition can occurs solely if the vapor is supersaturated 
(. S s > 0). It means there exist a A P driving the atoms toward the substrate. The eq. C.l 
summarizes this condition: 

AG V = ln(l + S s ), (C.l) 

Va 

where A G v is the change in the free energy per unit volume; Kb is the Boltzmann constant; 
V a the atomic volume; and S s = Pv p Ps is the vapor supersaturation with Pv-Ps being the 
difference between pressures of the vapor and above the solid, respectively. 

The molecules (or atoms) impinge onto the substrate and produce, in a reasonable ap¬ 
proximation, spherical nucleus of mean radius r. The change AG for the substrate-film sys¬ 
tem is due to three contributions: ( 1 ) the first associated with all bonds to create the nucleus 
(~ ? ,:i AC,.); (2) a superficial film/vapor term owing the new nucleus (~ r 2 "//,,); (3) and the 
contribution emerging from the difference between the energy at the old substrate/vapor surface 
(~ r 2 j sv ) and at the new substrate/film interface (~ r 2r ) s f). It reads: 

AG = air 3 AG v + a 2 r 2 j fv + a 3 (r 2 j fs - r 2 -f sv ), (C.2) 

where a\ = |(2 — 3 cosO + cos 3 6)\ a 2 = 27r(l — cos0)\ a 3 = irsen 2 9', and 6 is the wetting 
angle. 

In the fig. C.3 there is a plot of AG as function of r. At an instable equilibrium point, 
= 0 3 , there is a critical nucleus size, r*, from which island growth becomes energetically 
favorable. In other words, at least an amount of AG* = AG(r*) energy must be supplied to the 
system for triggering spontaneous growth, where the nucleus size grows indefinitely. 

An important question raises up: How islands reach the stable size , since for r < r* 
they tend to shrink? One must to take into account that theses processes are stochastic and 
thermal fluctuations ( KbT ss AG*) are responsible to yield a density N* of stable islands per 

3 And ^4^ < 0. 
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unit time. Assuming a Boltzmann statistic, then N* ~ exp[— AG*/KbT]. Since AG* must 
be function of T and of the molecular flux F [atoms /cm 2 s], one can get the feeling about how 
these parameters affect nucleation processes. 



Figure C.3: Qualitative change of the free energy as function of the mean nucleus 
size during a gas-to-solid phase transition. 


C.2.1 Nucleation dependence on substrate temperature 

It is reasonable that molecular flux is directly related to S such that one can replace 1 + S s by 
F/F e in the eq. C.l, with F e being the evaporation rate from the film at equilibrium. Inserting 
this in AG and extracting r* from the instable equilibrium point, one can show (considering 
d'y/dT « 0 N/mK [1]): 


(8r*/dT) lF ~ [a 2lfv + o 3 (7/a - lsv)\/T 2 . (C.3) 

At this point it is clear that the behavior of the critical size, as function of T, depends on 
the relationship between the surface tensions. Indeed, this relationship is also used to describe 
growth modes of films (see next section). For instance, taking the simplest situation, where 
7/s ~ 7 sv, one has (dr*/dT )| F > 0. Flence, locally, mean critical size of the islands becomes 
larger as temperature increases (at least for a range of T, because the result r\_/ 1/T 2 ). By using 
the same procedure one can also find: 

d(AG*)/dT) lp > 0, (C.4) 

likewise r*, the energy barrier AG* increases as T increases. This result implies that the nu¬ 
cleation rate is reduced exponentially in accordance with the Boltzmann statistic. Qualitatively, 
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theses results have already been found in the classical work on supercooled liquid tin and wa¬ 
ter drops, backing to the 1950’s [177]. The crystallizations of synthetic granite and granodior- 
ite [178], as well as the famous growth of Si on Si(001) [179] agree with those results. However, 
there are also examples going in the opposite direction. S. O. Ferreira et al. [8,9] have explored 
direct three-dimensional island formation in the growth of CdTe on Si(l 11) substrates for show¬ 
ing that the Quantum-Dots (QD’s) density increases with T in the range of 200 to 300 °C, whilst 
the mean QD’s size decreases (fig. C.4). 



Figure C.4: Atomic force microscopy images of CdTe Quantum-Dots grown by HWE 
on Si(111) substrates at (a) 200 °C, (b) 250 °C, and (c) 300 °C. (d) Dot density (■) and 
diameter (•) as function of T. Courtesy of Prof. S. O. Ferreira. These results can be 
found in the Ref. [8,9]. 


C.2.2 Nucleation dependence on molecular flux 

Once one has r* = r* (T. F ), it is straightforward showing the dependence on F at T constant. 
Considering the same simplest situation, one reaches to: 

(3r*/dF) w < 0 and [d(AG*)/dF\ ]T < 0. (C.5) 

Eq. C.5 tells us that nuclei size are smaller as the molecular flux increases, but at the 
same time, the nucleation rate is bigger. This result agree with the intuitive feeling that when 
more particles arrive onto substrate per unit time there is no time enough for the particles find 
an already stable island where they can pile up. Consequently, it is likely they form dimers 
and other instable small islands. These results match, for instance, with those reported for the 
growth of InAs/GaAs QD’s on GaAs(001) substrates by MBE [182]. At T = 490 °C and 
F = 0.016 ML s~\ the QD’s density is about 0.2 x 10 11 cm 2 and their diameter is near from 
200 A. However, as F is increased to 0.094 ML s -1 , at fixed temperature, the QD’s density also 
increases at least five times, whereas the QD’s diameter goes to about 150 A. 

Studies converging in this vein (and related to them) can be found in the cap. 9 of the 
Ref. [159] as well as in the already cited Refs. [8,9, 177-179, 182] and Refs, therein. 


Dot diameter ( nm) 








Appendix C - Concepts for Surfaces in Equilibrium and for the Thin Film Growth 


81 


C.3 Growth and structure of films 

C.3.1 Growth Modes 

The mechanical equilibrium condition (eq. C. 6 ), predicted by the drop liquid model, provide a 
reasonable way to explain three basic growth modes for solids. 

7 sv = 7sf + Ifv cos 9. (C. 6 ) 

When the film “wets” the whole substrate surface, nucleating two-dimensional islands 
on the interface, then 9 = 0 and y sv > 7 s f + 77 ,. This is the called Frank-van derMerwe (FM) 
or layer-by-layer growth mode, often observed in metal deposited on metal and homoepitaxial 
growth [1], In this mode, ideally, the subsequent layer starts growing just after the complete 
formation of the precedent layer. 

Otherwise, when there is not the “wetting layer” and three-dimensional islands are nucle¬ 
ated directly on the substrate, the Volmer-Weber (VW) growth mode takes place. In this situation 
parts of the substrate remain exposed and the energy is minimized reducing the film/vapor inter¬ 
face, since 7 f v ^$> y sv + 7 s j. CdTe grown on Si(l 11) (fig. C.4), for instance, follows this kind 
of growth [8,9]. 

Lastly, the Stranski-Krastanov (SK) growth mode is understood as a transition between 
FM to VW mode [1], At the first 3-5 monolayers of deposition, y sv dominates growth mak¬ 
ing wetting. After that, direct three-dimensional islands start to nucleate on the subsequent 
layers. Any energetic disturb (defects, accumulated strain, etc.) can be the origin of this 
transition [1-3]. A typical example of SK growth mode occurs in Ge films deposited on 
Si(001), where a metastable phase containing small clusters precedes the “macroscopic” three- 
dimensional growth. See the fig. C.5 [183], 

It is very important mention that the growth mode depends not only on the materials 
involved in the process, but mainly on the growth parameters, including whether impurities are 
present. The growth of Ag on Ag(lll) in the range of temperature from 175 to 575 K, for 
instance, changes gradually from three-dimensional (VW) to step-flow growth 4 . However, if 
submonolayer deposits of Sb are present, so the growth behavior becomes layer-by-layer (FM) 
over the entire range of 225K to 375 K [184]. 

Explained in the section 2.3.3. 
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Figure C.5: Ge films grown on Si(001) through SK mode, (a) Scanning Tunneling 
Microscopy (STM) 2500 x 2500 A image of “hut” Ge clusters on Ge layers grown in 
the FM growth mode. The clusters are < 1000 A long and 20-40 A high. One believes 
that this is a signature preceding macroscopic three-dimensional nucleations. (b) SEM 
image of the subsequent growth dominated by the macroscopic three-dimensional Ge 
islands. These clusters are ~ 250 A high and much larger than the hut clusters, which 
can not be viewed in this scale. Images extracted and edited from [183]. 


C.3.2 Commensurability and Polycrystallinity 

During a nucleation, under very special conditions, the film can “copy” the substrate crystalline 
structure [1], In these cases, one says that the growth is commensurate or epitaxial [3], Several 
semiconductor devices depends on epitaxy. The most commons are optoelectronic ones, as 
LEDs and lasers based on GalnN materials, and those used in high-speed microelectronic digital 
(wireless communications), composed of Ge.,.Si i _. r [1]. Typically, the first step in epitaxial 
films is the homo-epitaxial Si on Si growth by CVD. The reason is that the epilayer is free from 
defects, purer than the wafer and can be doped independently of it. 

Epitaxy depends on the mismatch between the lattice parameters (|ay — a s \/a s ) of the 
film (a f) and of the substrate (a s ), including their thermal dependence via the dilatation coeffi¬ 
cient [150, 155, 156], It occurs whether the relative mismatch is below ^ 10 - 15% [155,156], 
Nevertheless, exotic cases (called hard epitaxy) as the hetero-epitaxial CdTe layers grown on 
Si( 111) also exist [153], in spite of almost 19% and 48% mismatch of the lattice parameter and 
thermal coefficient, respectively. Larger mismatch implies larger accumulated strain (elastical 
energy) which can be released as defects or can support a transition from FM to VW growth 
mode. 

It is worth mention that the effect of heteroepitaxy can be included in the drop liquid 
model replacing the term (~ A G v ) by [~ (AG„ + AG S )] in the eq. C.2, where A G s inserts 
the extra energy accumulated in the strain form. AG* increases leading the nucleation rate to 
decline. Other effects as impurities, defects, chemical reactions also can be inserted into the 
model in order to get the felling about their consequences on nucleations [1]. 

Unlike epitaxial single-crystals, polycrystalline films are composed by an collection of 
grains, each of them having its own crystallographic orientation. These environments lead to a 
complex competition between neighboring grains which are droved by the surface tensions, once 
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7 = 7 (hkl) [1-3]. Generally, a preferential direction of growth (A), called texture [1, 150, 185], 
appears and evolves as deposition proceeds. GrainsA grow faster than the others and, eventually, 
dominate the surface. Texture affects mechanical and electrical features of films, including their 
elastic modulus, yield strength, magnetic permeability, etching rate, diffusion rate, and others 
[1]. Standard procedure quantifying texture is based on X-ray diffraction techniques [150,185]. 


C.4 Kinetic Phenomena and Superficial Structures 

Crystal (surface) growth is a far-from-equilibrium phenomenon. During a growth from vapor 
phase, atoms impinge onto the substrate (interface) breaking and/or forming local bonds, chang¬ 
ing the morphology and hampering the system to minimize its free energy [159]. As the sub¬ 
strate is heated, the adosrbed atoms, or molecules, (called adatoms ) might diffuse on the surface 
performing Brownian motions. Commonly, these adatoms stick at positions which maximize 
then - number of coordination and deposition occurs [1]. Thermal fluctuation may also induces 
desorption, i.e, adatoms leave the surface, returning to the vapor [22], Deposition, desorption 
and aggregation mediated by diffusion are the basic mechanisms ruling a general growth. Of 
course, specific situations excludes one, maybe two of these processes, but we are considering 
the wide case. These last two, and other cited below, are thermally activated, namely, they occur 
with rates given by Arrhenius’ laws (eq. C.7) [1-3,22]: 

r = t 0 exp [~E t /K b T\, (C.7) 

where r is a rate of a particular event, to a constant, and E r is the energy associated with the t 
process. 



A crystal surface is not a perfect plane. Indeed, it presents large pieces ( terraces ) sep- 
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V 



Figure C.7: Schematic of the lattice potential. 


arated by steps of atomic height (or multiples of an atomic heigth) - see and follow the fig. 
C.6. Even theses steps are not perfectly straight, rather they contain straight parts separated by 
kinks [3]. Particles arriving at the substrate (a) have an amount of energy relative to the va¬ 
por phase [2]. They use this energy for breaking chemical bonds already established at surface 
and/or for becoming adatoms (b). Adatoms diffuse with a rate since that, for each hop, 

an adatom must overcome the lattice potential existing between neighboring positions (see fig. 
C.7). This potential depends on the nature of compounds involved in the growth and on the 
substrate face. 

Preferential sites for sticking are the edge of steps and kinks, once they have a larger 
number of dangling bonds [3]. However, adatoms reach these positions only if the temperature 
is sufficiently high for promoting a length diffusion ( Xdiff ) comparable to the terrace separation 
length ( xterr ) [2]. Step-flow growth is achieved when Xdiff 3> xt err [3]. Interestingly, if an 
adatom diffuses on a terrace, but towards another one (c), it tends to come back instead to 
jump ojf the terrace. This is known as Ehrlich-Schowoebel ejfect (ES) [159, 180]. In the hg. 
C.7 the effect of the ES barrier is added. In cases where the ES barrier dominates the growth 
(see Pt grown on Pt(lll), hg. 102 in the ref. [159]) structures looking like wedding cakes are 
formed [159]. Among several ways to define an expression for the diffusion energy that includes 
the ES barrier, one usually employed is: 


^diff — Eo T nE n + E es , 


(C.8) 


where n is the number of coordination; E n is the energy per unit bond and E es is the associated 
energy to the ES bander. Typically, E n ~ 0.1 eV and E es ~ 0.01 eV [22]. However, it is worth 
keep in mind that these values may change appreciably depending on the system considered. 
Additionally, the eq. C.8 is solely a simplified model which does not catch all source of energies 
involved in real processes. 

Since an adatom establishes bonds, deposition happens. When an adatom meets another 
atom, they can form a dimer (d). Dimers also diffuse, albeit it is (much)less probable due to 
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the enhance of Ediff. Sometimes, dimers grow and become a two- or three-dimensional island 
(e), which tend to capture adatoms (f) (if r > r*). At the same time, some atoms break off 
then bonds and leave the surface (g). So, desoiption occurs. Desorption is an effect competing 
with deposition. Due to this effect, the molecular flux migth be different from the deposition 
rate. However, in usual MBE (and related technique) conditions, desorption is a neglegible 
process [1-3,22]. For Ga on GaAs(l 11) over the range 860K to 960 K, the desorption energy is 
Edes ~ 2.5 eV [22], Assuming l/rodes = 10~ 14 s and T = 850 K in the inverse of eq. C.7 5 , 
one obtains that only one Ga atom is desorbed from the surface in each 2 s. Comparing this 
desorption rate with the deposition rate used in MBE (~ 1 0 16-18 atoms cm -2 s _1 ), it eases to 
see that desorption phenomenon is inoperative [ 22 ]. 


^The inverse of a given rate is the characteristic time of this process. 
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Wsat 

Saturation value for the roughness 

w{L , t ) 

Global roughness 

W 

Thermodynamical work 

x dif f 

Characteristic length for diffusion 

X dif f 

Length of a terrace 

Z 

Dynamic exponent 


Greek Symbols 

Description 

a 

Roughness exponent 

OL 1 

Geometrical exponent 


Spectral roughness exponent 

7 

Surface energy or Surface tension 

7 fv 

Surface tension at the film-vapor interface 

7p 

Non-universal parameter in the KPZ scaling theory context 

7 sf 

Surface tension at the substrate-film interface 

'7 sv 

Surface tension at the substrate-vapor interface 

r 

Non-universal parameter in the KPZ scaling theory context 

r / 

Gamma function 

r (i,t) 

Slope-Slope spatial covariance 

e 

Tilting angle 

c 

Characteristic length in the surface 

Cp 

Non-universal parameter in the KPZ scaling theory context 

V 

Noise 

Vp 

Non-universal parameter in the KPZ scaling theory context 

®hkl 

Wetting angle 

Ohkl 

Angle of reflection for the (hkl) peak in a XRD measurement 

0 

Functional ofx, h and t 

K 

Kappa exponent 

A 

Excess of velocity of an interface 

^CuKa 

Wavelenght of the radiation in a K a transition for a Cu atom 

A 

Preferential direction of growth - Texture direction 
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I^V 

V 

t 

gi 


? 

T 

TOdes 

T diff 

$ 

X 

(x> 
(x 2 )c 
Ipk 

0 


Chemical potential 
Chemical potential in a vapor phase 
Strength of the “surface tension” 
Correlation length 
Parallel correlation length 
Grand potential 
Standard deviation ofh 
Standard deviation of x 
Scale factor 

Rate of a particular kinetic effect 

Rate of desorption 

Rate of diffusion 

Universal scaling function 

Random variable 

Mean of chi 

Second cumulant of x 

Poly gamma function of kth-order 

Universal scaling function 

Aspect ratio of a grain/mound 
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